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ABSTRACT 


Current two-dimensional preliminary design codes use structured programming, 
which is rigid and does not allow the user to vary parameters easily. This study uses 
object-oriented programming to allow the user to vary all selectable parameters in a 
familiar Windows operating environment. The programmed design is based on the 
assumptions of axial and free-vortex flow between blade rows, simple radial equilibrium, 
and a thermally and calorically perfect gas. The program allows a fan or core stage 
design and uses an open architecture to facilitate upgrades and extensions. 

Using the Naval Postgraduate School’s (NPS) transonic compressor design as 
input, the preliminary design code output was compared to the detailed throughflow 
design of the transonic compressor. The results agreed reasonably well with detailed 
throughflow design. With some minor improvements this code can easily be used to 


develop a preliminary design that can be optimized to the user’s requirements. 
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I. INTRODUCTION 


The design of a new axial compressor involves a sequence of steps, progressing 
through a sequence of computational programs of increasing complexity and 
sophistication. The first, or ‘preliminary design’ step, can be a one-dimensional ‘mean- 
line calculation’, or a two-dimensional calculation of a preliminary flow path and 
selection of the blading. The latter is the minimum required if the overall task is the 
preliminary design of an aircraft gas turbine engine. It is also what is required to provide 
the input for detailed throughflow design codes [Ref. 1], which, in turn, generate inputs to 


codes, which compose the blade geometry [Ref. 2] for manufacturing. 


The current preliminary compressor design code used in aircraft engine design 
courses [Ref. 3], was developed progressively using different versions of Hewlett- 
Packard (HP) BASIC [Ref. 4]. This highly structured programming language is rigid and 
does not allow the user to vary one, or several, parameters easily in order to change or 
optimize a design. Use of the programs requires the installation of HP BASIC for 


Windows. 


The purpose of the present study was to develop a preliminary compressor design 


code that would satisfy the following conditions: 
e Be simple to use. 
e Allow all selectable parameters to be changed. 
e Require only the Windows operating system. 


e Use an open architecture to allow upgrades (e.g., different loss models) or 


additions (e.g., turbine design) 


To meet these conditions, Microsoft’s Visual Basic 6.0 [Ref. 5] was selected and 


used to develop the source code. 
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Il. PROGRAM OVERVIEW 


A. ASSUMPTIONS 


The following assumptions were made in the development of the programmed 


equations and design of the code: 


e Design is for an axial flow compressor (with radial movement of the mean 


line). 
e Simple radial equilibrium is assumed from hub to tip. 
e Free vortex flow is assumed between blade rows. 


e A conceptual engine design study [using Ref. 6 for example] will generate the 


inputs for the code. 


B. INPUTS AND OUTPUTS 


The stage is shown schematically in Fig. 1, and the inputs required by the code, 


are shown in Table 1. 






Mass Flow 
Total Pressure 
Total Temperature 


Stage Pressure Ratio 
Stage Temperature Ratio 
Stage Efficiency 


—=> 


ae in: 


o 


[i CENTERLINE 


Figure 1. Throughflow Schematic 


Inlet (1) Stator 3) 





Table 1. User-Provided Inputs 


At the inlet station, it is required to specify the gas (through the gas constant and 
ratio of specific heats; a thermally and calorically perfect gas is assumed), the mass flow 
rate, and the stagnation conditions. Then, in order to accommodate the usual fan or core 
design constraints, four parameters must be specified for the rotor; namely, rotational 
speed, axial Mach number, relative Mach number, and flow angle at the tip. Only three 


of these four are independent. Off-line calculation is required to determine the other one. 


The design selections for the rotor and stator are listed in columns two and three 
of Table 1. Note that the blade loading throughout the stage is determined by the 
selection of the diffusion factor at only one location for one blade. Blade aspect ratio and 
thickness variations are choices, which are determined by structural considerations. 


Structural constraints are not included in the program. 


The parameters, which are calculated and output by the code, are shown in 


Table 2. 


tet) Retort ator 
Comp. | Hub | Mean | Tip Hub | Mean |_ Tip Hub | Mean | Tip 
(h) (m) (m) (m) 
Xih Xim Xit X2m X3m 
Xz1m XzZ2m XzZ3m 
Xuim Xvu2m Xyu3m 
Xwim Xw2m Xw3m 
Xoim Xe2u Xo3m 
Yim 12m 13m 
M lm Mom M3m 
Mzim Mz2m Mz3m 
Mwim Mw2m Mw3m 
B 1m Bom Bam 
Oim O2m O3m 
Yom Y3 m 
Yw2m Yw3m 
Input D3m 
Input Input 
Tst2m Tst3m 
Le 


2m 3m 


0*om 0*3m 


é mn ) ae 
Thr Tht3 

Ao A3 
Ty/Tr1 T3/Te1 
P2/P2 P3/P 3 
T/T T3/Tu 





Table 2. Outputs 
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Il. ALGORITHMS 


A. PROGRAM SEQUENCE 


The overall program sequence is shown in Fig. 2. After the program is started a 
“splash” screen is first shown. A splash screen is an introductory screen (similar to an 
“about” screen) which states the name, owner(s), and version of the program. Next, an 
interactive screen appears, giving the user a choice of either going through a compressor 
design or a turbine design. The scope of this paper covers only a compressor design; 
therefore the turbine option has been disabled. When the Compressor Design button is 
pressed the input screen will appear. The design sequence is shown in Fig. 3. The 
required inputs are shown in Table 1. Once the inputs are typed in the user presses the 
OK button and the main screen appears. The main screen is laid out in a tab format with 


five tabs. The tabs are as follows: 
e Inlet Conditions (1) 
e Rotor Calculations (2) 
e Stator Calculations (3) 
e Stage Performance 
e Blade Geometry 


This allows the variables and code to be grouped so the user can better understand 
what is being displayed. It also allows the programmer to develop the code with an open 
architecture for easy updates. At this point the user can still modify the inputs before 


executing the design calculations. 


START 
PROGRAM 









CODE NOT 
AVAILABLE 






SCREEN 
SHOWING 
OPTIONS FOR 
TURBINE OR 
COMPRESSOR 
DESIGN 









COMPRESSOR 





SCREEN SHOWING 
REQUIRED INPUTS 
(TABLE 1) 











USER INPUTS 
CONCEPTUAL 
DESIGN OUTPUT 





MAIN SCREEN 
DISPLAYED 









USER MODIFIES 
INPUTS (IF 
NECESSARY) 





Figure 2. 





CODE NOT 
TURBINE AVAILABLE ) NEXT 
STAGE 









CALCULATE BASIC 
THROUGHFLOW OR 

CALCULATE INPUTS FOR 
NEXT STAGE? 







APPLY INPUTS 


CALCULATE 
STAGE 


THROUGHFLOW 
& GEOMETRY 









INPUT SCREEN FOR 
THE WHOLE NUMBER 
OF ROTOR AND 
STATOR BLADES 








MODIFY 
INPUTS 


OUTPUT 
DISPLAYED 
ON MAIN 
SCREEN 













EXIT, MODIFY 
INPUTS OR NEXT 
STAGE? 






EXIT 


END PROGRAM 


Main Program Flowchart 


INPUT 
PARAMETERS 


CALCULATE 
INLET 
CONDITIONS 


CALCULATE VELOCITY 
DIAGRAMS ACROSS 
ROTOR MEAN-LINE 


CALCULATE VELOCITY 
DIAGRAMS ACROSS 
STATOR MEAN-LINE 


CALCULATE ROTOR & 
STATOR LOSSES 


CALCULATE ROTOR & 
STATOR PERFORMANCE 
& ANNULUS 


CALCULATE HUB & TIP 
VELOCITY DIAGRAMS 
FOR ROTOR & STATOR 


CALCULATE INITIAL 
BLADE GEOMETRY 





Figure 3. 








REQUEST USER TO 
INPUT WHOLE 
NUMBER OF ROTOR 
& STATOR BLADES 


CALCULATE REVISED 
BLADE GEOMETRY 


RECALCULATE 
LOSSES 


NOTE: SHOCK, 
SECONDARY FLOW AND 
TIP CLEARANCE LOSSES 
ARE INITIALLY SET TO 0 


DO THE LOSSES 
CONVERGE? 


NO 


YES 


CALCULATE CAMBER, 
INCEDENCE, & 
DEVIATION ANGLE 


INTERPOLATION 
SUBROUTINE 
[REF 10-11] 


CALCULATE 
COMPRESSOR 
EFFICIENCY 





DISPLAY 
OUTPUT 


Design Flowchart 
° 


The user has the option to either apply the inputs to the design calculations or 
calculate the next stage. A multistage approach is also outside the scope of this paper so 
this option has been disabled. After the user presses the Apply button the program 
calculates the basic throughflow velocity diagrams, losses, performance and geometry. 
During the throughflow calculations an input screen is displayed for the user to input the 
whole number of blades needed (Z) for the rotor and stator. The main screen is displayed 
again with the outputs filled in the text boxes. The outputs are shown in Table 2. With 
the initial throughflow calculations complete the user may end the program of modify the 


inputs until the desired results are achieved. 


B. DESIGN EQUATIONS 


The equations programmed in the source code are listed in Appendix A. The 
equations are grouped the same as the tabs on the main screen. The equations for the 
velocity diagrams, profile loss, secondary flow loss, tip clearance loss, and stage 
performance are from Shreeve [Ref. 7 and 8]. The shock loss equation is from Koch and 


Smith [Ref. 9]. 


The incidence, deviation and camber angles (which relate the flow angles to the 
blade geometry, as shown in Fig. 4) were derived from NASA SP-36 [Ref. 10]. Sixth 
degree polynomial curvefits were used to approximate the data in Figs. 137, 138, 142, 
161, 162, 168, 172, 178, 179, and 180 of Ref. 10. An interpolation routine [Ref. 11], was 


programmed to solve for unknowns within the curve-fits. 
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Figure 4. Blade Geometry Schematic 


The performance of the stage is calculated following the thermodynamic process 
shown in Fig. 5. The conditions on the mean line are taken as being representative of the 
stage; however, the loss coefficients include contributions due to secondary flow and tip- 


clearance. 
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IV. PROGRAM STRUCTURE 


A. OBJECTS 


Forms (a.k.a. screens) are used to interact with the user. Forms are ‘objects’ in 
Visual Basic and Visual C++ programming. Other common objects include text boxes, 
option boxes and tabs. Objects have properties (this is how the object looks to the user) 
which can be affected during design or run time (program execution). The objects allow 
the user to modify the inputs throughout the program and either executes the desired 
modification or keeps the original inputs. The objects (or “screen captures”) are shown 


in Appendix B.1. 


B. VARIABLES 


Variables are used in programming as placeholders of data where they are used 
throughout the code to execute statements. Two-dimensional arrays are used extensively 
throughout the code since they allow a parameter that changes both throughflow and 
spanwise to be stored in one variable. For example, Xw (relative velocity) varies from 
inlet (1) to stator (3) and from hub (1) to tip (3). If we store this variable separately we 
would need to track nine variables. Instead, we simply track one variable Xw(3,3) where 
the first number is the throughflow number and the second number is the spanwise 
number. Appendix B.2 shows all the variables used in the program. This includes all the 


input and output parameters as well as others that are not seen by the user. 


C. MODULES 


Code is used to state what needs to be executed. The code syntax is based on 
BASIC, which translates “pseudo” language to machine language. Code is written for 
objects to tell the object what to do when acted upon by the user. For this program there 
is very little code for objects. Most of the code is a sequence of design equations which 
do not cause changes in any objects and require minimal interaction with the user during 
execution. Modules are used in this situation. Modules take subroutines or functions that 


do not affect an object and keep them in a separate file for easier reading during 
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programming or debugging. The code, under the form ‘frmCompressor’, is the design 
code sequence shown in Figure 2. Appendix B.3 shows all the source code for the 


program. 
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V. RESULTS AND DISCUSSION 


Sanger’s transonic compressor design [Ref. 13] was used to test the ability of the 
code to approximate a known axial stage design. A set of hand calculations was also 
carried out independently of the code in order to both validate the coding, and to 
document the test case. The design input flow conditions, and parameters derived from 
the final geometry of the Sanger design, were used as inputs to both the preliminary 
design code and hand calculations. Appendix C.1 and Appendix C.2 document, in detail, 


the results of the hand calculations and the preliminary design code, respectively. 


The results from hand calculations and from the code were compared to the 
output of the streamline curvature code applied to the Sanger design, which is given in 
Appendix E of Sanger [Ref. 13]. The comparisons are shown in Appendix D. It can be 
seen from the comparison charts in Appendix D that the hand calculations agreed fully 
with the preliminary design code calculations. Also, for most parameters, the hand and 
preliminary design code calculations agreed with the streamline curvature code outputs. 
In Figure D1, the calculated annulus geometry agrees well at the inlet and then begins to 
deviate somewhat through the stage. This is because the streamline curvature code takes 
into account blockage, whereas the preliminary design code and hand calculations, do 
not. This is easy to correct. The stage performance however, shown in Fig. D2, is 
predicted very successfully by the code. In Fig. D3 and D4, the differences in solidity and 
blade height are also, indirectly, the result of omitting blockage from the calculation of 
annulus area. Velocity diagram details are compared in Figs. D5 to D8. It is clear that the 


preliminary design code reproduces the final design values to very acceptable accuracy. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


From the comparisons given in Appendix D, it can be seen that the preliminary 
design code does very well in developing the velocity diagrams and the initial blading 
geometry necessary for a detailed throughflow design and final geometry calculation 


[Ref. 3 and 4]. 


Improvements can be made in order to have the design more detailed as well as 


improve the code’s usability. They are as follows: 


e Add the ability to do a multistage design (use the output of the previous stage 


as an input to the new stage). 
e Add different input screens for specific design cases (e.g., fan or core). 
e Add stress limits (hoop and centrifugal) for fan design. 
e Add blockage and bleed to the throughflow calculations. 
e Draw scale velocity diagrams. 
e Draw blades based on code output (1.e., built in geometry package). 
e Add parametric analysis for a range of values (e.g., inlet flow angle). 
e Add turbine stage design (in parallel). 
e Add the ability to open and save data. 
e Add the ability to print the user’s results. 
e¢ Compile the code into a stand-alone executable. 


Incorporation of these improvements will make the code a preliminary 
turbomachine design software package that can be used as inputs to detailed design 


packages, as well as providing a needed teaching tool for aircraft engine design. 
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APPENDIX A. DESIGN EQUATIONS 


A.1 shows the throughflow equations and A.2 shows the interpolation equations 
for the incidence, camber and deviation angles. 


A.l. THROUGHFLOW EQUATIONS 
1. Inlet Conditions 


Given: m,P,,T, Mwi @,R, Y 
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3. Stator Across the Mean Line 
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4. Rotor Performance 
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5. Stator Performance 
Initially set @..,@. =0 


On subsequent iterations @,,,,.. and @, are 
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6. Rotor & Stator at the Hub and Tip 


The below equations are generic forms of the specific equation. The brackets [] replace the rotor 
(2) and stator (3) subscripts. Also, the parentheses () replace the hub (1) and tip (3) subscripts. 
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The below equations are generic forms of the specific equation. The brackets [] replace the rotor 
(2) and stator (3) subscripts. 
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A.2. INTERPOLATION EQUATIONS 


Starting with a general quadratic equation 
D'(R)=A,-R*+B,-R+C 


solve for A;, B;, and C,; using known points i-1, 1, and i+1. 
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APPENDIX B. PROGRAM SOURCE DATA 
B.1. SCREEN SNAPSHOTS 


1 sy) 
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Figure B2. 
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Figure B4. 





Compressor Design Screen: Inlet Conditions ran 
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Figure B6. Compressor Design Screen: Stator Calculations (3) 
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Figure B8. Compressor Design Screen: Blade Geometry | 
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Blade Number Input Screen — 


Figure B9. 





About Screen 


Figure B10. 
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B.2. KEY VARIABLES 


Error! Not a valid link. 


Table B1. Velocity Diagram Variables 


Error! Not a valid link. 


Table B2. Stage Performance Variables 


Error! Not a valid link. 
Table B3. Blade Geometry Variables 
B.3. SOURCE CODE 


Lines with an apostrophe (') in the front is a comment line in the code and is not 
an executed statement. 


1. Splash Screen Code 


Option Explicit 
Private Sub Form_Load() 
IbIVersion.Caption = "Version " & App.Major & "." & App.Minor & "." & App.Revision 
Ib]ProductName.Caption = App.Title 
Ibl1Company.Caption = App.CompanyName 
End Sub 
Private Sub tmrSplash_Timer() 
Unload Me 
frmCompTurb.Show 


End Sub 
2. Compressor/Turbine Selection Screen Code 


Private Sub cmdComp_Click() 
Unload Me 
frmInitial.Show 


End Sub 
3: Initial Input Screen Code 


Option Explicit 
Option Base | 
Private Sub OKButton_Click() 


frmInitial. Hide 
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frmCompressor.Show 
' Transfer initial values to the main form. 


End Sub 
4. Compressor Design Screen Code 


Option Explicit 

Option Base | 

Public i As Integer 

‘Integer to use as an array counter from inlet (i = 1) to rotor (i=2) to stator (i=3) 


‘It should be noted that although some arrays have no inlet values (e.g. diffusion and deg. of 
reation) 


‘we still set the array value to 3 versus 2. This is to keep consistent with the numbering. 
‘arrays that don't have the i=1 value assigned are null and irrelevent since they are not used 
‘in any calculations. 

Public j As Integer 

‘ integer to use as an array counter from hub (j = 1) to mean (j = 2) to tip (j = 3) 

Tt should be noted that although hub, mn and tip have been set up a constants 


'j will be used to calculate the minimum loss incidence angle and the camber angle 


Public ksh As Single 
‘correction factor for shape 
Public slopegraph As Boolean 
‘checks which graph should be used to calculate the slope factor 
Const testvar As Single = 0.0001 
"exit criteria for main loop 
Public loopcount As Integer 
‘loop counter to track the number of iterations 
Dim tempvar(3) As Double 
holder for previous loss value for comparison 
Public CompEff As Double 
‘compressor efficiency 
Private Sub cmdApply_Click() 
Sequence 
Update VelDiag 
UpdatePerf 
UpdateGeo 
End Sub 
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Private Sub Form_InitializeQ 
frmCompressor.txtmdot.Text = frmInitial.txtmdot.Text 
frmCompressor.txtPtl.Text = frmInitial.txtPt1.Text 
frmCompressor.txtTt1.Text = frmInitial.txtTtl.Text 
frmCompressor.txtR.Text = frmInitial.txtR.Text 
frmCompressor.txtMzt(1).Text = frmInitial.txtMz1t.Text 
frmCompressor.txtMwt(1).Text = frmInitial.txtMw 1t.Text 
frmCompressor.txtspeed.Text = frmInitial.txtomega.Text 
frmCompressor.txtgamma.Text = frmInitial.txtgamma.Text 
frmCompressor.txtalphat(1).Text = frmInitial.txtAlphalt.Text 
frmCompressor.txtDm(2).Text = frmInitial.txtDRm.Text 
frmCompressor.txtsigmam(2).Text = frmInitial.txtsigmaRm.Text 
frmCompressor.txtsigmam(3).Text = frmInitial.txtsigmaSm.Text 
frmCompressor.txtR21.Text = frmInitial.txtR21.Text 
frmCompressor.txtR32.Text = frmInitial.txtR32.Text 
frmCompressor.txtphi21.Text = frmInitial.txtphi21.Text 
frmCompressor.txtphi32.Text = frmInitial.txtphi32.Text 
frmCompressor.txtA31.Text = frmInitial.txtA31.Text 
frmCompressor.txtDelta(2).Text = frmInitial.txtDeltar.Text 
frmCompressor.txtDelta(3).Text = frmInitial.txtDeltas.Text 
frmCompressor.txtAR(2).Text = frmInitial.txtARr.Text 
frmCompressor.txtAR(3).Text = frmInitial.txtARs.Text 
frmCompressor.txttch(2).Text = frmInitial.txttcrh.Text 
frmCompressor.txttch(3).Text = frmInitial.txttcsh.Text 
frmCompressor.txttcm(2).Text = frmInitial.txttcrm.Text 
frmCompressor.txttcm(3).Text = frmInitial.txttcs m.Text 
frmCompressor.txttct(2).Text = frmInitial.txttcrt.Text 
frmCompressor.txttct(3).Text = frmInitial.txttcst.Text 
frmCompressor.txtfsigma(2).Text = frmInitial.txtfsigmar.Text 
frmCompressor.txtfsigma(3).Text = frmInitial.txtfsigmas.Text 
"Initialize Inlet Condition Variables. 
alpha(1, 3) = frmCompressor.txtalphat(1).Text 
gamma = frmCompressor.txtgamma.Text 
Mz(1, 3) = frmCompressor.txtMzt(1).Text 
mdot = frmCompressor.txtmdot.Text 
Mw(1, 3) = frmCompressor.txtMwt(1).Text 


omega = frmCompressor.txtspeed.Text 
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Pt(1) = frmCompressor.txtPt!.Text 

Rbar = frmCompressor.txtR.Text 

Tt(1) = frmCompressor.txtTtl.Text 

‘Initialize Rotor Assumption Variables 

D(2, 2) = frmCompressor.txtDm(2).Text 
sigma(2, 2) = frmCompressor.txtsigmam(2).Text 
pitch(2) = frmCompressor.txtR21.Text 

phi21 = frmCompressor.txtphi21.Text 

fsigma(2) = frmCompressor.txtfsigma(2).Text 
‘Initialize Stator Assumption Variables 

A31 = frmCompressor.txtA3 1.Text 

sigma(3, 2) = frmCompressor.txtsigmam(3).Text 
pitch(3) = frmCompressor.txtR32.Text 

phi32 = frmCompre ssor.txtphi32.Text 

fsigma(3) = frmCompressor.txtfsigma(3).Text 
‘Initialize Stage Performance variables 

delta(2) = frmCompressor.txtDelta(2).Text 
delta(3) = frmCompressor.txtDelta(3).Text 


"initialize blade geometry variables 


For i= 2 To3 
For j = 1 To3 
ARG) = frmCompressor.txtAR(i).Text 
If j = 1 Then 
tc(i, J) = frmCompressor.txttch(i).Text 
Elself j = 2 Then 
tc(i, j) = frmCompressor.txttcm(i).Text 
ElseIf j = 3 Then 
tc(i, j) = frmCompressor.txttct(i).Text 
End If 
Next j 
Next i 
End Sub 
Private Sub UpdateVelDiag() 
For i=1To3 
For j = 1 To3 
If} = 1 Then 


txtBetah(i).Text = beta(i, j) 
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End If 


txtAlphah(i).Text = alpha(i, j) 
txtXh(i).Text = x(i, j) 
txtXUh(i). Text = Xu(i, j) 
txtXzh(i).Text = Xz(i, j) 
txtXthetah(i).Text = Xtheta(i, j) 
txtXwh(i).Text = Xw(i, j) 
txtrh(i).Text = RG, j) 
txtMh(i).Text = MG, j) 
txtMwh(i).Text = Mwii, j) 
txtMzh(i).Text = Mz(i, j) 
If i= 2 Then 
txtrsth(i).Text = rst(i, j) 
txtDh(i).Text = DG, j) 
txtYh(i).Text = Y(i, j) 
txtY wh(i).Text = Yw(i, j) 
Elself i = 3 Then 
txtrsth(i).Text = rst(i, j) 
txtDh(i).Text = DG, j) 
txtYh(i).Text = Y(i, j) 
txtY wh(i).Text = Ywi(i, j) 
End If 


If } = 2 Then 


txtbetam(i).Text = beta(i, j) 
txtalpham(i).Text = alpha(i, j) 
txtXm(i).Text = x(i, j) 
txtXUm(i).Text = Xu(i, j) 
txtXzm(i).Text = Xz(i, j) 
txtXthetam(i).Tex = Xtheta(i, j) 
txtXwm(i).Text = Xw(, j) 
txtrm(i).Text = R(i, j) 
txtMm(i).Text = M(i, j) 
txtMwm(i).Text = Mw, j) 
txtMzm(i).Text = Mz(i, j) 
If i= 2 Then 

txtrstm(i).Text = rst(i, j) 

txtYm(i).Text = Y(i, j) 
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txtYwm(i).Text = Yw@, j) 


Elself i = 3 Then 


End If 
End If 
If } = 3 Then 


txt rstm(i).Text = rst(i, j) 
txtYm(i).Text = Y(i, j) 
txtYwm(i).Text = YwG, j) 
txtDm(i) = DG, j) 


txtbetat(i).Text = beta(i, j) 


txtXt(i).Text = x(i, j) 


txtXUt(i).Text = Xui, j) 
txtXzt(i).Text = Xz(i, j) 


txtXthetat(i).Text = Xtheta(, j) 


txtXwt(i).Text = Xw(i, j) 


txtrt(i). Text = R(i, j) 
txtMt(i). Text = M(i, j) 
If i= 2 Then 


txtalphat(i).Text = alpha(i, j) 
txtMwt(i).Text = MwG, j) 
txtMzt(i).Text = Mz(i, j) 
txtrstt(i).Text = rst(i, j) 
txtDt(i).Text = D(i, j) 

txt Yt(i).Text = Y(i, j) 

txtY wt(i).Text = Yw(i, j) 


Elself i = 3 Then 


End If 
End If 


txtalphat(i).Text = alpha(i, j) 
txtMwt(i).Text = MwG, j) 
txtMzt(i).Text = Mz(i, j) 
txtrstt(i).Text = rst(i, j) 
txtDt(i).Text = D(i, j) 

txt Yt(i).Text = Y(i, j) 
txtYwt(i).Text = Yw(i, j) 


txtA(i). Text = A(i) 


txtrht(i). Text = rht(i) 
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Next j 
Next i 
txttau.Text = tau 
End Sub 
Private Sub UpdatePerf() 
For i=1To3 
txtTTt1(i).Text = TTt1(i) 
txtPPt1(i).Text = PPt1(i) 


Next i 

For i=1To2 
txtPEPt1(i).Text = PEPt1 (i) 

Next i 

For 1 = 2 To 3 
txtDelta(i).Text = delta(i) 
txtomegasftc(1).Text = omegasftc(1) 
txtomegap(1).Text = omegap(i) 
txtomegasl(i).Text = omegas(i) 
txtomega(i). Text = omegat(i) 
txtAA1(i).Text = AAI(i) 
txtPhim(i).Text = Capphi(i, mn) 
txtPtPt1(i).Text = PtPt1(i) 

Next i 


txtTETt1(1).Text = TETt1(1) 
txtTR1ITt1.Text = Tr1Ttl 
txtPR1Pt1.Text = Pri Ptl 
txtTt3Ttl Text = TtTt1(3) 
txtPI.Text = PtPt1(3) 
txtBTau.Text = tau 
txteff.Text = CompEff 
End Sub 
Private Sub UpdateGeo() 
For 1 = 2 To 3 
txtH().Text = hii) 
txtC(i).Text = Ci) 
txtZ(i).Text = Z(i) 
txtZrev(i).Text = intZ(i) 
txtARrev(i).Text = ARrev(i) 
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txtCrev(i).Text = Crev(i) 
txtS().Text = Sq) 
txtcamberh(i) = camber(i, 1) 
txtcamberm(i) = camber(i, 2) 
txtcambert(i) = camber(i, 3) 
txtincidenceh(i) = icor(i, 1) 
txtincidencem(i) = icor(i, 2) 
txtincidencet(i) = icor(i, 3) 
txtdevh() = dref(i, 1) 
txtdevm(i) = dref(i, 2) 
txtdevt(i) = dref(i, 3) 


Next i 


Private Sub Form_Load() 


‘Initialize Inlet Condition Variables. 

alpha(1, 3) = frmCompressor.txtalphat(1).Text 
gamma = frmCompressor.txtgamma.Text 
Mz(1, 3) = frmCompressor.txtMzt(1).Text 
mdot = frmCompressor.txtmdot.Text 

Mw(1, 3) = frmCompressor.txtMwt(1).Text 
omega = frmCompressor.txtspeed.Text 

Pt(1) = frmCompressor.txtPt!.Text 

Rbar = frmCompressor.txtR.Text 

Tt(1) = frmCompressor.txtTtl .Text 

‘Initialize Rotor Assumption Variables 

D(2, 2) = frmCompressor.txtDm(2).Text 
sigma(2, 2) = frmCompressor.txtsigmam(2).Text 
pitch(2) = frmCompressor.txtR21.Text 

phi21 = frmCompressor.txtphi21.Text 
fsigma(2) = frmCompressor.txtfsigma(2).Text 
‘Initialize Stator Assumption Variables 

A31 = frmCompressor.txtA3 1.Text 

sigma(3, 2) = frmCompressor.txtsigmam(3).Text 
pitch(3) = frmCompressor.txtR32.Text 

phi32 = frmCompressor.txtphi32.Text 
fsigma(3) = frmCompressor.txtfsigma(3).Text 


‘[nitialize Stage Performance variables 


4] 


delta(2) = frmCompressor.txtDelta(2).Text 
delta(3) = frmCompressor.txtDelta(3).Text 


"initialize blade geometry variables 


For 1 = 2 To 3 
For j =1To3 
AR(i) = frmCompressor.txtAR(i).Text 
If j = 1 Then 
tc(i, J) = frmCompressor.txttch(i).Text 
Elself j = 2 Then 
tc(i, j) = frmCo mpressor.txttcm(i).Text 
Elself j = 3 Then 
tc(i, J) = frmCompressor.txttct(i).Text 
End If 
Next j 
Next i 


End Sub 

Private Sub mnuAbout_Click() 
frmAbout.Show 

End Sub 

Private Sub mnuContents_Click() 
MsgBox "This function is not available yet.", vbInformation, "Not Implemented" 

End Sub 

Private Sub mnuExit_Click() 
End 

End Sub 

Private Sub mnuOpen_Click() 
digFile.ShowOpen 

End Sub 

Private Sub mnuPrint_Click() 
digFile.ShowPrinter 
Printer.Copies = dlgFile.Copies 
Printer.Orientation = dlgFile.Orientation 
Printer.Print A (i) 
Printer.EndDoc 

End Sub 

Private Sub mnuSave_Click() 
digFile.ShowSave 


42 


End Sub 
Private Sub mnuSearch_Click() 
MsgBox "This function is not available yet.", vbInformation, "Not Implemented" 
End Sub 
Public Sub Sequence() 
i=l 
InletCond 
HubCalc i 
If optNACA = True Then 
ksh = 1 
slopegraph = True 
End If 
If optC = True Then 
ksh = 1.1 
slopegraph = False 
End If 
If optDCA = True Then 
ksh = 0.7 
slopegraph = False 


End If 

For 1 = 2 To 3 
MeanCalc i 

Next i 


loopcount = 1 
Do 
For i=2 To3 
If loopcount = | Then 
omegat(i) = 0 
omegasftc(i) = 0 
omegas(i) = 0 
End If 
tempvar(i) = omegat(i) 
If i= 2 Then 
RotorPerf i, loopcount 
ElseIf i = 3 Then 
StatorPerf i, loopcount 


End If 
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HubTipCalc i 
BladeGeo i, loopcount 
Next i 
loopcount = loopcount + 1 


Loop Until Abs(omegat(2) - tempvar(2)) < testvar And Abs(omegat(3) - tempvar(3)) < 
testvar 


For i=2 To3 
For j = 1 To3 
Incidence i, j, slopegraph, ksh 
Next j 
Next i 
CompEff = ((PtPt1(3) * glg(gamma)) - 1) / (tau - 1) 
End Sub 
Private Sub txtA31_KeyPress(KeyAscii As Integer) 
Select Case KeyAscii 
Case vbKey0 To vbKey9 
Case vbDecimal, 46 
Case vbKeyBack, vbKeyClear, voKeyTab, vbKeyUp, vbKeyDown 
Case 45 
If Len(txtA31.Text) <> 0 Then 
KeyAscii = 0 ' ignore keystroke 
Beep 
End If 
Case Else 
KeyAscii = 0 
Beep 
End Select 
End Sub 
Private Sub txtalphat_KeyPress(Index As Integer, KeyAscii As Integer) 
Select Case KeyAscii 
Case vbKey0 To vbKey9 
Case vbDecimal, 46 
Case vbKeyBack, vbKeyClear, voKeyTab, vbKeyUp, vbKeyDown 
Case 45 
If Len(txtalphat(1).Text) <> 0 Then 
KeyAscii = 0 ' ignore keystroke 
Beep 


44 


End If 


Case Else 
KeyAscii = 0 
Beep 
End Select 


End Sub 
Private Sub txtAR_KeyPress(Index As Integer, KeyAscii As Integer) 
Select Case KeyAscii 

Case vbKey0 To vbKey9 

Case vbDecimal, 46 

Case vbKeyBack, vbKeyClear, voKeyTab, vbKeyUp, vbKeyDown 

Case 45 

If Len(txtAR(2).Text) <> 0 Then 
KeyAscii = 0 ' ignore keystroke 


Beep 

ElseIf Len(txtAR(3).Text) <> 0 Then 
KeyAscii = 0 
Beep 

End If 

Case Else 
KeyAscii = 0 
Beep 
End Select 
End Sub 


Private Sub txtDelta_KeyPress(Index As Integer, KeyAscii As Integer) 
Select Case KeyAscii 
Case vbKey0 To vbKey9 
Case vbDecimal, 46 
Case vbKeyBack, vbKeyClear, voKeyTab, vbKeyUp, vbKeyDown 
Case 45 
If Len(txtDelta(2).Text) <> 0 Then 
KeyAscii = 0 ' ignore keystroke 


Beep 

ElseIf Len(txtDelta(3).Text) <> 0 Then 
KeyAscii = 0 
Beep 


End If 
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Case Else 
KeyAscii = 0 
Beep 
End Select 
End Sub 
Private Sub txtDm_KeyPress(Index As Integer, KeyAscii As Integer) 
Select Case KeyAscii 
Case vbKey0 To vbKey9 
Case vbDecimal, 46 
Case vbKeyBack, vbKeyClear, voKeyTab, vbKeyUp, vbKeyDown 
Case 45 
If Len(txtDm(2).Text) <> 0 Then 
KeyAscii = 0 ' ignore keystroke 


Beep 
End If 
Case Else 
KeyAscii = 0 
Beep 
End Select 
End Sub 


Private Sub txtgamma_KeyPress(KeyAscii As Integer) 
Select Case KeyAscii 
Case vbKey0 To vbKey9 
Case vbDecimal, 46 
Case vbKeyBack, vbKeyClear, vbKeyTab, vbKeyUp, vbKeyDown 
Case 45 
If Len(txtgamma.Text) <> 0 Then 


KeyAscii = 0 ' ignore keystroke 


Beep 
End If 
Case Else 
KeyAscii = 0 
Beep 
End Select 
End Sub 


Private Sub txtmdot_KeyPress(KeyAscii As Integer) 
Select Case KeyAscii 
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Case vbKey0 To vbKey9 
Case vbDecimal, 46 
Case vbKeyBack, vbKeyClear, voKeyTab, vbKeyUp, vbKeyDown 
Case 45 
If Len(txtmdot.Text) <> 0 Then 


KeyAscii = 0 ' ignore keystroke 


Beep 
End If 
Case Else 
KeyAscii = 0 
Beep 
End Select 
End Sub 


Private Sub txtMwt_KeyPress(Index As Integer, KeyAscii As Integer) 
Select Case KeyAscii 
Case vbKey0 To vbKey9 
Case vbDecimal, 46 
Case vbKeyBack, vbKeyClear, voKeyTab, vbKeyUp, vbKeyDown 
Case 45 
If Len(txtMwt(1).Text) <> 0 Then 


KeyAscii = 0 ' ignore keystroke 


Beep 
End If 
Case Else 
KeyAscii = 0 
Beep 
End Select 
End Sub 


Private Sub txtMzt_KeyPress(Index As Integer, KeyAscii As Integer) 
Select Case KeyAscii 
Case vbKey0 To vbKey9 
Case vbDecimal, 46 
Case vbKeyBack, vbKeyClear, voKeyTab, vbKeyUp, vbKeyDown 
Case 45 
If Len(txtMzt(1).Text) <> 0 Then 
KeyAscii = 0 ' ignore keystroke 
Beep 
47 


End If 


Case Else 
KeyAscii = 0 
Beep 
End Select 


End Sub 
Private Sub txtphi21_KeyPress(KeyAscii As Integer) 
Select Case KeyAscii 

Case vbKey0 To vbKey9 

Case vbDecimal, 46 

Case vbKeyBack, vbKeyClear, voKeyTab, vbKeyUp, vbKeyDown 

Case 45 

If Len(txtphi21.Text) <> 0 Then 
KeyAscii = 0 ' ignore keystroke 


Beep 
End If 
Case Else 
KeyAscii = 0 
Beep 
End Select 
End Sub 


Private Sub txtphi32_KeyPress(KeyAscii As Integer) 
Select Case KeyAscii 
Case vbKey0 To vbKey9 
Case vbDecimal, 46 
Case vbKeyBack, vbKeyClear, voKeyTab, vbKeyUp, vbKeyDown 
Case 45 
If Len(txtphi32.Text) <> 0 Then 
KeyAscii = 0 ' ignore keystroke 


Beep 
End If 
Case Else 
KeyAscii = 0 
Beep 
End Select 
End Sub 


Private Sub txtPt!_KeyPress(KeyAscii As Integer) 
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Select Case KeyAscii 
Case vbKey0 To vbKey9 
Case vbDecimal, 46 
Case vbKeyBack, vbKeyClear, voKeyTab, vbKeyUp, vbKeyDown 
Case 45 
If Len(txtPt1.Text) <> 0 Then 
KeyAscii = 0 ' ignore keystroke 


Beep 
End If 
Case Else 
KeyAscii = 0 
Beep 
End Select 
End Sub 


Private Sub txtR_KeyPress(KeyAscii As Integer) 
Select Case KeyAscii 
Case vbKey0 To vbKey9 
Case vbDecimal, 46 
Case vbKeyBack, vbKeyClear, voKeyTab, vbKeyUp, vbKeyDown 
Case 45 
If Len(txtR.Text) <> 0 Then 
KeyAscii = 0 ' ignore keystroke 


Beep 
End If 
Case Else 
KeyAscii = 0 
Beep 
End Select 
End Sub 


Private Sub txtR21_KeyPress(KeyAscii As Integer) 
Select Case KeyAscii 
Case vbKey0 To vbKey9 
Case vbDecimal, 46 
Case vbKeyBack, vbKeyClear, voKeyTab, vbKeyUp, vbKeyDown 
Case 45 
If Len(txtR21.Text) <> 0 Then 
KeyAscii = 0 ' ignore keystroke 
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Beep 


End If 
Case Else 
KeyAscii = 0 
Beep 
End Select 
End Sub 


Private Sub txtR32_KeyPress(KeyAscii As Integer) 
Select Case KeyAscii 
Case vbKey0 To vbKey9 
Case vbDecimal, 46 
Case vbKeyBack, vbKeyClear, voKeyTab, vbKeyUp, vbKeyDown 
Case 45 
If Len(txtR32.Text) <> 0 Then 


KeyAscii = 0 ' ignore keystroke 


Beep 
End If 
Case Else 
KeyAscii = 0 
Beep 
End Select 
End Sub 


Private Sub txtsigmam_KeyPress(Index As Integer, KeyAscii As Integer) 
Select Case KeyAscii 
Case vbKey0 To vbKey9 
Case vbDecimal, 46 
Case vbKeyBack, vbKeyClear, voKeyTab, vbKeyUp, vbKeyDown 
Case 45 
If Len(t xtsigmam(2).Text) <> 0 Then 
KeyAscii = 0 ' ignore keystroke 
Beep 
ElseIf Len(txtsigmam(3).Text) <> 0 Then 


KeyAscii = 0 
Beep 
End If 
Case Else 
KeyAscii = 0 
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Beep 
End Select 
End Sub 
Private Sub txtspeed_KeyPress(KeyAscii As Integer) 
Select Case KeyAscii 
Case vbKey0 To vbKey9 
Case vbDecimal, 46 
Case vbKeyBack, vbKeyClear, voKeyTab, vbKeyUp, vbKeyDown 
Case 45 
If Len(txtspeed.Text) <> 0 Then 
KeyAscii = 0 ' ignore keystroke 


Beep 
End If 
Case Else 
KeyAscii = 0 
Beep 
End Select 
End Sub 


Private Sub txttch_KeyPress(Index As Integer, KeyAscii As Integer) 
Select Case KeyAscii 
Case vbKey0 To vbKey9 
Case vbDecimal, 46 
Case vbKeyBack, vbKeyClear, voKeyTab, vbKeyUp, vbKeyDown 
Case 45 
If Len(txttch(2).Text) <> 0 Then 
KeyAscii = 0 ' ignore keystroke 


Beep 

ElseIf Len(txt tch(3).Text) <> 0 Then 
KeyAscii = 0 
Beep 

End If 

Case Else 
KeyAscii = 0 
Beep 
End Select 
End Sub 


Private Sub txttcm_KeyPress(Index As Integer, KeyAscii As Integer) 
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Select Case KeyAscii 
Case vbKey0 To vbKey9 
Case vbDecimal, 46 
Case vbKeyBack, vbKeyClear, vbKeyTab, vbKeyUp, vbKeyDown 
Case 45 
If Len(txttem(2).Text) <> 0 Then 
KeyAscii = 0 ' ignore keystroke 
Beep 
ElseIf Len(txttcm(3).Text) <> 0 Then 


KeyAscii = 0 
Beep 

End If 

Case Else 
KeyAscii = 0 
Beep 
End Select 
End Sub 


Private Sub txttct_KeyPress(Index As Integer, KeyAscii As Integer) 
Select Case KeyAscii 
Case vbKey0 To vbKey9 
Case vbDecimal, 46 
Case vbKeyBack, vbKeyClear, voKeyTab, vbKeyUp, vbKeyDown 
Case 45 
If Len(txttct(2).Text) <> 0 Then 
KeyAscii = 0 ' ignore keystroke 
Beep 
ElseIf Len(txttct(3).Text) <> 0 Then 


KeyAscii = 0 
Beep 

End If 

Case Else 
KeyAscii = 0 
Beep 
End Select 
End Sub 


Private Sub txtTtl_KeyPress(KeyAscii As Integer) 
Select Case KeyAscii 
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Case vbKey0 To vbKey9 
Case vbDecimal, 46 
Case vbKeyBack, vbKeyClear, voKeyTab, vbKeyUp, vbKeyDown 
Case 45 
If Len(txtTt1.Text) <> 0 Then 


KeyAscii = 0 ' ignore keystroke 


Beep 
End If 
Case Else 
KeyAscii = 0 
Beep 
End Select 
End Sub 


5. Blade Number Input Screen Code 


Option Explicit 
Option Base | 
Private x As Variant 
Private Sub OKButton_Click() 
X = txtZ.Text 
If IsNumeric(x) = False Then 
MsgBox "Input is not a number!", vbExclamation, "Numeric Validation” 
Else 
frmZ.Hide 
End If 
End Sub 


6. About Screen Code 


Option Explicit 

"Reg Key Security Options... 

Const READ_CONTROL = &H20000 

Const KEY_QUERY_VALUE = &H1 

Const KEY_SET_VALUE = &H2 

Const KEY_CREATE_SUB_KEY = &H4 

Const KEY_ENUMERATE_SUB_KEYS = &H8 
Const KEY_NOTIFY = &H10 

Const KEY_CREATE_LINK = &H20 
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Const KEY_ALL_ACCESS = KEY_QUERY_VALUE + KEY_SET_VALUE + 
_ KEY_CREATE_SUB_KEY + KEY_ENUMERATE_SUB_KEYS + 
_ KEY_NOTIFY + KEY_CREATE_LINK + READ_CONTROL 

"Reg Key ROOT Types... 

Const HKEY_LOCAL_MACHINE = &H80000002 

Const ERROR_SUCCESS = 0 

Const REG_SZ= 1 ‘Unicode nul terminated string 

Const REG_DWORD = 4 "32-bit number 

Const gREGKEYSYSINFOLOC = "SOFTWARE\Microsoft\Shared Tools Location" 

Const gREGVALSYSINFOLOC = "MSINFO" 

Const gREGKEYSYSINFO = "SOFTWARE\Microsoft\Shared Tools \MSINFO" 

Const gREGVALSYSINFO = "PATH" 


Private Declare Function RegOpenKeyEx Lib "advapi32" Alias "RegOpenKeyExA" (ByVal hKey 
As Long, ByVal IpSubKey As String, ByVal ulOptions As Long, ByVal samDesired As 
Long, ByRef phkResult As Long) As Long 


Private Declare Function RegQueryValueEx Lib "advapi32" Alias "RegQueryValueExA" (ByVal 
hKey As Long, ByVal IpValueName As String, ByVal IpReserved As Long, ByRef 
IpType As Long, ByVal IpData As String, ByRef lpcbData As Long) As Long 


Private Declare Function RegCloseKey Lib "advapi32" (ByVal hKey As Long) As Long 
Private Sub cmdSysInfo_ClickQ) 
Call StartSysInfo 
End Sub 
Private Sub cmdOK_Click() 
Unload Me 
End Sub 
Private Sub Form_Load() 
Me.Caption = "About " & App.Title 
IblVersion.Caption = "Version " & App.Major & "." & App.Minor & "." & App.Revision 
IbITitle.Caption = App.Title 
IblDescription.Caption = App.FileDescription 
End Sub 


Public Sub StartSysInfo() 
On Error GoTo SysInfoErr 
Dim re As Long 
Dim SysInfoPath As String 
‘Try To Get System Info Program Path\Name From Registry... 
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If GetKey Value(HKEY_LOCAL_MACHINE, gREGKEYSYSINFO, 
gREGVALSYSINFO, SysInfoPath) Then 


‘Try To Get System Info Program Path Only From Registry... 


Elself GetKey Value(HKEY_LOCAL_MACHINE, gREGKEYSYSINFOLOC, 
gREGVALSYSINFOLOC, SysInfoPath) Then 


‘Validate Existance Of Known 32 Bit File Version 
If (Dir(SysInfoPath & "\MSINFO32.EXE") <> "") Then 
SysInfoPath = SysInfoPath & "\MSINFO32.EXE' 
Error - File Can Not Be Found... 
Else 
GoTo SysInfoErr 
End If 
‘Error - Registry Entry Can Not Be Found... 
Else 
GoTo SysInfoErr 
End If 
Call Shell(SysInfoPath, vbNormalFocus) 
Exit Sub 
SysInfoErr: 
MsgBox "System Information Is Unavailable At This Time", voOKOnly 
End Sub 


Public Function GetKeyValue(KeyRoot As Long, KeyName As String, SubKeyRef As String, 
ByRef KeyVal As String) As Boolean 


Dim i As Long "Loop Counter 

Dim re As Long "Return Code 

Dim hKey As Long ‘Handle To An Open Registry Key 

Dim hDepth As Long ' 

Dim KeyValType As Long ‘Data Type Of A Registry Key 

Dim tmpVal As String 'Tempory Storage For A Registry Key Value 


Dim KeyValSize As Long ' Size Of Registry Key Variable 


rc = RegOpenKeyEx(KeyRoot, KeyName, 0, KEY_ALL_ACCESS, hKey) 


"Open Registry Key 
If (rc <> ERROR_SUCCESS) Then GoTo GetKeyError ‘Handle Error... 
tmpVal = String$(1024, 0) ' Allocate Variable Space 
KeyValSize = 1024 "Mark Variable Size 
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rc = RegQueryValueEx(hKey, SubKeyRef, 0, _ 


KeyValType, tmpVal, KeyValSize) ‘ Get/Create Key Value 
If (rc <> ERROR_SUCCESS) Then GoTo GetKeyError ‘Handle Errors 
If (Asc(Mid(tmpVal, KeyValSize, 1)) = 0) Then "  Win95 Adds Null 
Terminated String... 
tmpVal = Left(tmpVal, KeyValSize - 1) ‘Null Found, Extract 
From String 
Else ‘WinNT Does NOT Null Terminate String... 
tmpVal = Left(tmpVal, KeyValSize) ‘Null Not Found, 
Extract String Only 


Select Case KeyValType ' Search Data Types... 

Case REG_SZ ' String Registry Key Data Type 
KeyVal = tmpVal Copy String Value 

Case REG_DWORD ‘Double Word Registry Key Data Type 


For i = Len(tmpVal) To | Step—1 ' Convert Each Bit 
KeyVal = KeyVal + Hex(Asc(Mid(tmpVal, i, 1))) ' Build Value Char. 


By Char. 

Next 

KeyVal = Format$("&h" + KeyVal) ‘Convert Double Word To String 
End Select 
GetKey Value = True "Return Success 
rc = RegCloseKey(hKey) Close Registry Key 
Exit Function " Exit 

GetKeyError: ‘Cleanup After An Error Has Occured... 

KeyVal ="" Set Return Val To Empty String 
GetKey Value = False Return Failure 
rc = RegCloseKey(hKey) Close Registry Key 


End Function 


Module mdlInletCond Code 


Option Explicit 
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Option Base | 

Public gamma As Double 

' Specific Heat Ratio (Cp/Cv) 
Public mdot As Double 
"Mass Flow 

Public Rbar As Double 

' Specific Gas Constant 
Public omega As Double 
‘wheel speed 

Public Cp As Double 

‘at constant pressure 

Const g As Long = 32.2 

' gravitational constant 
Public Capphi(3, 3) As Double 
"Flow Function 

Public Tt(3) As Double 
‘Total Temperature 

Public Pt(3) As Double 

' Total Pressure 

Public M(3, 3) As Double 
"Mach Number 

Public Mw(3, 3) As Double 
"Mach Relative to the Blade 





Public Mz(3, 3) As Double 





Mach of the Axial Component 
Public alpha(3, 3) As Double 

‘Inlet Flow angle 

Public beta(3, 3) As Double 

‘Inlet Flow angle 

Public x(3, 3) As Double 
‘Dimensionless Velocity 

Public Xz(3, 3) As Double 
‘Dimensionless Velocity of the Axial Component 
Public Xtheta(3, 3) As Double 
‘Dimensionaless Velocity along Theta 
Public Xu(3, 3) As Double 


‘Dimensionless Velocity of the Wheel Speed Component 
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Public Xw(3, 3) As Double 
‘Dimensionless Velocity Relative to the Blade 
Public rot(3) As Double 
"Density 
Public Vt(3) As Double 
' Total Velocity 
Public A(3) As Double 
‘Annulus Area 
Public R(3, 3) As Double 
' Radius 
Public rht(3) As Double 
Hub to Tip Ratio 
Public pitch(3) As Double 
"Mean Line Pitch (Rm2/Rm1) 
Public rst(3, 3) As Double 
‘Degree of Reaction 
Public i As Integer 
Public Const hub As Integer = | 
‘Integer used when making hub calculations (hub = 1) 
Public Const mn As Integer = 2 
‘Integer used when making mean line calculations (mn = 2) 
Public Const tip As Integer = 3 
‘Integer used when making tip calculations (tip = 3) 
Public Sub InletCond() 
Calculation of Inlet Conditions at the tip 
i=l 
beta(i, tip) = Arccos(Mz(i, tip) / Mwd, tip)) 
MA, tip) = Mz(i, tip) / Cos(DegToRad(alphati, tip))) 
x(i, tip) = Sqr((((gamma - 1) / 2) * MG, tip) * 2) / (1 + ((gamma - 1) / 2) * MG, tip) * 2))) 
Xz(i, tip) = x(i, tip) * Cos(DegToRad(alpha(i, tip))) 
Xtheta(i, tip) = Xz(i, tip) * Tan(DegToRad(alpha(i, tip))) 
Xu(i, tip) = Xtheta(i, tip) + Xz(i, tip) * Tan(DegToRad(beta(i, tip))) 
Xw(i, tip) = Xwfunc(Xz(i, tip), beta(i, tip)) 
rot(i) = Pt(i) / (12 * Rbar * Tt(i)) 
Cp = ggl(gamma) * Rbar 
Vt(i) = Sqr(2 * Cp * g * Tt(i)) * 12 
Capphi(i, tip) = xi, tip) * (1 - xG, tip) * 2) * (1 / (gamma - 1)) 
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A(i) = (mdot / (rot(i) * VtG))) * (1 / (CapphiGi, tip) * Cos(DegToRad(alpha(i, tip))))) 
RG, tip) = (XuG, tip) * Vt(i)) / omega 
RG, hub) = Sqr(RG, tip) * 2 - AG) / (22/7)) 
rht(i) = rhtfunc(RG, tip), RG, hub)) 
"Mean Calculations 
RG, mn) = rmfunc(RG, tip), RG, hub)) 
Xtheta(i, mn) = Xthetafunc(R(i, tip), RG, mn), Xtheta(i, tip)) 
Xz(i, mn) = Xz(i, tip) 
alpha(i, mn) = alphafunc(Xtheta(i, mn), Xz(i, mn)) 
Xu(i, mn) = Xufunc(R(, tip), RG, mn), Xu, tip)) 
beta(i, mn) = betafunc(Xu(i, mn), Xtheta(i, mn), Xz(i, mn)) 
x(i, mn) = Xfunc(Xz(i, mn), alpha(i, mn)) 
Mi, mn) = Mach(gamma, x(i, mn)) 
Xw(i, mn) = Xwfunc(Xz(i, mn), beta(i, mn)) 
Mw(i, mn) = Mach(gamma, Xw(i, mn)) 
Mz(i, mn) = Machz(M@, mn), alpha(i, mn)) 
End Sub 


8. Module mdlHubCalc Code 


Option Explicit 
Option Base | 
Public sigma(3, 3) As Double 
' Solidity 
Public Sub HubCalc(i As Integer) 
Hub Calculations 
Xz(i, hub) = Xz(i, mn) ' Radial Equilibrium 
Xtheta(i, hub) = Xthetafunc(R(, mn), RG, hub), Xtheta(i, mn)) 
alpha(i, hub) = alphafunc(Xtheta(i, hub), Xz(i, hub)) 
x(i, hub) = Xfunc(Xz(i, hub), alpha(i, hub)) 
Xu(i, hub) = Xufunc(RG, mn), RG, hub), Xu, mn)) ' for constant wheel speed 
beta(i, hub) = betafunc(Xu(i, hub), Xtheta(i, hub), Xz(i, hub)) 
Xw(i, hub) = Xwfunc(Xz(i, hub), beta(i, hub)) 
Md, hub) = Mach(gamma, x(i, hub)) 
Mw, hub) = Mach(gamma, Xw(i, hub)) 
Mz(i, hub) = Machz(M(i, hub), alpha(i, hub)) 
End Sub 
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9, Module mdlMeanCalc Code 


Option Explicit 
Option Base | 
Public Y(3, 3) As Double 
‘Dimensionless Local Velocity 
Public Yw(3, 3) As Double 
‘Dimensionless Local Velocity Relative to the Blade 
Public Yz(3, 3) As Double 
‘Dimensionless Local Velocity along the Axial Component 
Public A31 As Double 
Ratio of alpha. Unity for a repeating stage. 
Public D(3, 3) As Double 
' Diffusion Factor 
Public tau As Double 
Total temperature ratio 
Public phi21 As Double 
‘Change in Axial Velocity across the Rotor. 
Public phi32 As Double 
‘Change in Axial Velocity across the Stator. 
Public phi2m As Double 
‘Change in Axial Velocity at the Rotor exit at the Mean Line. 
Public philm As Double 
‘ratio of axial velocity to rotation velocity 
Public Sub MeanCalc(i As Integer) 
‘Mean Calculations for the Rotor or Stator 
If i= 2 Then 
Dim temp As Double 


beta(i, mm) = Arcsin(SinB2(DiffA(sigma(i, mn), pitch(i), DG, mn), phi21, beta(i 
- 1, mn)), DiffB(sigma(i, mn), pitch(i)), pitch(i))) 


Xu(i, mn) = pitch(i) * XuGi - 1, mn) 

philm = Xz@i - 1, mn) / XuG@i - 1, mn) 

phi2m = phi21 * philm * (1 / pitch(i)) 

Xz(i, mn) = phi2m* Xu, mn) 

Xtheta(i, mn) = Xu(i, mn) - Xz(i, mn) * Tan(DegToRad(beta(i, mn))) 

alpha(i, mn) = alphafunc(Xtheta(i, mn), Xz(i, mn)) 

tau = taufunc(Xu(i, mn), Xtheta(i, mn), Xu(i- 1, mn), Xtheta(i - 1, mn)) 
Else 
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alpha(i, mn) = A31 * alpha(i - 2, mn) 


DQG, mn) = Diffusion(phi32, alpha’ - 1, mn), alpha(i, mn), pitch(i), sigma(i, 
mn)) 


Xz(i, mn) = phi32 * Xz(i- 1, mn) 
Xtheta(i, mn) = Xz(i, mn) * Tan(DegToRad(alpha(i, mn))) 
Xu(i, mn) = Xu(i- 1, mn) * pitch) 
beta(i, mn) = betafunc(Xu(i, mn), Xtheta(i, mn), Xz(i, mn)) 
End If 
rst(i, mn) = DofReaction(Xtheta(i - 1, mn), Xu(i- 1, mn), Xtheta(i, mn), Xu(i, mn)) 
x(i, mn) = Xfunc(Xz(i, mn), alpha(i, mn)) 
Xw(i, mn) = Xwfunc(Xz(i, mn), beta(i, mn)) 
Y(i, mn) = Yfunc(x(i, mn), tau) 
Ywi, mn) = Yfunc(Xw(i, mn), tau) 
M(i, mn) = Mach(gamma, Y(i, mn)) 
Mw(i, mn) = Mach(gamma, Yw(i, mn)) 
Mz(i, mn) = Machz(M(i, mn), alpha(i, mn)) 
End Sub 


10. Module mdlStgPerformance 


Option Explicit 

Option Base 1 

Public omegasftc(3) As Double 
‘Secondary flow and tip clearance loss 
Public omegap(3) As Double 

' Profile Loss 

Public omegas(3) As Double 

' Shock Loss 

Public omegat(3) As Double 

' Total Loss (secondary flow + tip clearance + profile) 
Public TTt1(3) As Double 

' Static to Total Temperature Ratio 
Public PPt1(3) As Double 

' Static to Total Pressure Ratio 

Public Tr1Tt! As Double 

' Total Relative Temperature Ratio 
Public Pri Pt! As Double 


' Total Relative Pressure Ratio 
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Public TETt1(3) As Double 

Public PEPt1(3) As Double 

Public PtPt1(3) As Double 

Total Pressure Ratio 

Public TtTt1(3) As Double 

‘Total Temperature Ratio 

Public AA1(3) As Double 

' Area constriction ratio 

Public delta(3) As Double 

‘Tip Gap 

Public Sub RotorPerf(i As Integer, loopcount As Integer) 
Dim ShLoss As Double 


omegap(i) = 2 * sigma, mn) * ((Cos(DegToRad(betai - 1, mmn))) *% 2) / 
(Cos(DegToRad(beta(i, mn))) * 3)) * (0.005 + (0.16 * (DG, mn) * 4))) 


If loopcount > | Then 
If Mwd - 1, mn) > 1 Then 


ShLoss = ((1 - ShockLoss(Mw(i - 1, mn), gamma)) / (1 - (1 + ((gamma 
-1)/2) * Mw@- 1, mn) * 2) * -1 * ggl(gamma)))) 


End If 
omegas(1) = ShLoss 
omegasftc(i) = SFTC(betaGi - 1, mn), beta, mn), sigma(i, mn), h(i), delta(i), 


S@)) 
Else 
omegas(i) = 0 
omegasftc(i) = 0 
End If 


omegat(i) = omegap(i) + omegasftc(i) + omegas(1) 
If loopcount = | Then 
TTtl@ - 1)=1 -(xG- 1, mn) % 2) 
PPtl(@- 1)=TTtla - 1) * ggl(gamma) 
Tr1Ttl = TTtl@ - 1) + (Xw@i - 1, mn) * 2) 
PriPtl = Tr1Ttl * gg1(gamma) 
TETtl@ - 1) = Tr1Ttl + (XuG, mn) * 2) - (XuGi- 1, mn) * 2) 
PEPtl(i- 1) = TETt!1@ - 1) * ggl(gamma) 
End If 
PEPt1 (i) = PEPt1(i - 1) - (omegat(i) * (Pri Pt1 - PPtlq - 1))) 
TtTtl(@) = tau 
TETt1(i) = TETt1( - 1) 
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PtPt1 (i) = PEPt1(i) * ((TtTt1G) / TETt1(i)) * gg1(gamma)) ' Rotor Compression ratio 
TTtl() = TtTt1(@ - (x(i, mn) * 2) 

PPt1(i) = PtPtl(i) * ((TTt1(i) / tau) * gg1(gamma)) 

CapphiG, mn) = Y(i, mn) * ((1 - (YG, mn) * 2)) * g1(gamma)) 


AAI(i) = ((Capphidi - 1, tip) * Cos(DegToRad(alphaGi - 1, tip)))) / (CapphiG, mn) * 
Cos(DegToRad(alpha(i, mn))))) * ((Sqr(TtTt1(i))) / PtPt1()) 


A(i) = AAI(i) * AG - 1) 
RG, mn) = RG - 1, mn) * pitch(i) 
h(i) = AG) /(2 * (22/7) * RG, mn)) 
rht(i) = rhtfunc2(h@), RG, mn)) 
RG, tip) = (2/ (1 + rht(i))) * RG, mn) 
RG, hub) = rht(i) * RG, tip) 

End Sub 


Public Sub StatorPerf(i As Integer, loopcount As Integer) 
Dim ShLoss As Double 


omegap(i) = 2 * sigma(i, mn) * ((Cos(DegToRad(alphaGi - 1, mn))) *% 2) / 
(Cos(DegToRad(alpha(i, mn))) * 3)) * (0.005 + (0.16 * (DG, mn) * 4))) 


If loopcount > | Then 
If MG.- 1, mn) > 1 Then 


ShLoss = ((1 - ShockLoss(M(i - 1, mn), gamma)) / (1 - (1 + ((gamma - 
1)/2) * MG - 1, mn) * 2)“ (-1 * ggl(gamma)))) 


End If 
omegas(i) = ShLoss 
omegasftc(i) = SFTC(alphati - 1, mn), alpha(i, mn), sigma(i, mn), h(a), delta(i), 


S@) 
Else 
omegas(i) = 0 
omegasftc(i) = 0 
End If 


omegat(i) = omegap(i) + omegasftc(i) + omegas(1) 
PtPt1(i) = PtPt1(i - 1) - omegat(i) * (PtPt1(i - 1) - PPtl@ - 1)) 
TtTt1(i) = tau 
TTtl(i) = TtTt1(i) - (xG@, mn) * 2) 
PPt1(i) = ((TTt1@) / TtTt1@)) * ggl(gamma)) * PtPt1() 
Capphi(i, mn) = Yi, mn) * ((1 - (Yi, mn) * 2)) * g1(gamma)) 
AAI(i) = ((Capphii - 2, tip) * Cos(DegToRad(alpha@ - 2, tip)))) / (CapphiG, mn) * 
Cos(DegToRad(alpha(i, mn))))) * ((Sqr(TtTt1(i))) / PtPt1()) 
A(@i) = AAI(i) * AG - 2) 
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RG, mn) = RG - 1, mn) * pitch(i) 
h(i) = AG) / (2 * (22/7) * RG, mn)) 
rht(i) = rhtfunc2(h(@), RG, mn)) 
RG, tip) = (2/ (1 + rht(i))) * RG, mn) 
RG, hub) = rht(i) * RG, tip) 

End Sub 


11. Module mdlHubTipCalc Code 


Option Explicit 

Option Base | 

Dim 1 As Integer 

Public fsigma(3) As Double 


‘Linear function of sigma 
Public Sub HubTipCalc(i As Integer) 
‘Hub & Tip Calculations for the Rotor or Stator 
For 1 = 1 To 3 Step 2 
Xz(i, 1) = Xz(i, mn) 
Xtheta(i, 1) = (RG, mn) * Xtheta(i, mn)) / RG, 1) ' Radial Equilibrium 
Xu(i, 1) = Xufunc(RG, mn), RG, 1), XuG, mn)) ‘for constant wheel speed 
alpha(i, 1) = alphafunc(Xtheta(i, 1), Xz(i, 1)) 
beta(i, 1) = betafunc(XuG, 1), Xtheta(i, 1), Xz, 1)) 
x(i, 1) = Xfunc(Xzii, 1), alphaq, 1)) 
Xwi(i, 1) = Xwfunc(Xzii, 1), beta(i, 1) 
Y(i, J) = Yfunc(x(i, 1), tau) 
Ywid, 1) = Yfunc(Xwii, 1), tau) 
Ma, 1) = Mach(gamma, Y(i, 1)) 
Mwii, 1) = Mach(gamma, Ywdi, 1)) 
Mad, 1) = Machz(Mdi, 1), alphati, 1)) 
If 1= 1 Then 
sigma(i, 1) = sigma(i, mn) * (2 - fsigma(i)) 
Elself 1 = 3 Then 
sigma(i, 1) = sigma(i, mn) * fsigma(i) 
End If 
If i= 2 Then 
DG, 1) = Diffusion(phi21, beta(i - 1, 1), beta, 1), pitch(i), sigmati, 1)) 
Else 
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Da, 1) = Diffusion(phi32, alphadi - 1, 1), alpha(i, 1), pitch(i), sigmaci, 1)) 
End If 
rst(i, 1) = DofReaction(Xtheta(i - 1, 1), Xu - 1, 1), XthetaG, 1), Xu, 1) 
Next | 
End Sub 


12. Module mdlBladeGeometry Code 


Option Explicit 

Option Base | 

Public AR(3) As Double 

' Aspect Ratio 

Public C(3) As Double 

’ Chord 

Public Crev(3) As Double 

‘revised chord length 

Public ARrev(3) As Double 

‘revised aspect ratio 

Public Z(3) As Double 

"Number of Blades calculated 

Public intZ(3) As Integer 

"Number of Blades chosen by the user. 

Public h(3) As Double 

Blade height 

Public S(3) As Double 

blade spacing 

Public Sub BladeGeo(i As Integer, loopcount As Integer) 
h(i) = RG, tip) - RG, hub) ' blade height 
Ci) = hG) / ARG) ' Calculated chord length 
Z(i) = 2 * (22/7) * RG, mn) * sigma(i, mn) / C4) ' Calculated no. of blades 
If loopcount = | Then 

If i= 2 Then 


frmZ.1b1Z.Caption = "The calculated no. of blades in the rotor (Zr) is " 
& Z(i) 


Elself i = 3 Then 


frmZ.1b1Z.Caption = "The calculated no. of blades in the stator (Zs) is " 
& Z(i) 


End If 
frmZ.Show vbModal 
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intZ(i) = frmZ.txtZ.Text 
End If 
ARrev(i) = h(i) / (2 * (22/7) * RG, mn) * sigma(i, mn)) / intZ(i)) ' Revised aspect ratio 
frmZ.txtZ.Text = '"" 
Crev(i) = hG) / ARrev(i) ' Revised chord length 
SQ) = Crev(i) / sig ma(i, mn) ' Blade spacing 
End Sub 


13. Module mdlIncidence Code 


Option Explicit 
Option Base | 


Public 10(6, 9) As Double, slopen(7, 9) As Double, ikit(6) As Double, i0delta0(6, 9) As Double, 
slopemn(7, 9) As Double, Kit(7, 7) As Double, d0(6, 9) As Double, slopem(7, 9) As 
Double, dkit(6) As Double 


'X values for 6th order curve fit polynomial. 
Public 1010(3, 3) As Double 
'‘Zero-camber incidence angle 
Public d010(3, 3) As Double 
"Zero-camber deviation angle 
Public N(3, 3) As Double 
‘incidence angle slope factor n 
Public ikt(3, 3) As Double 
‘incidence angle correction factor for thickness 
Public dkt(3, 3) As Double 
‘ deviation angle correction factor for thickness 
Public i0ref(3, 3) As Double 
‘reference incidence angle 
Public dOref(3, 3) As Double 
‘reference deviation angle 
Public i2d As Double 
‘2 dimentional incidence angle 
Public icor(3, 3) As Double 
‘corrected 2 dimentional incidence angle 
Public dref(3, 3) As Double 
‘ deviation angle 
Public camber(3, 3) As Double 
"camber angle 
Public 10d0(3, 3) As Double 
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‘ variation of 10-d0 

Public ktbar(3, 3) As Double 

‘camber angle correction factor for thickness 

Public onemn(3, 3) As Double 

‘camber angle slope factor 1-m+n 

Public dm(3, 3) As Double 

‘ deviation angle slope factor m 

Public tc(3, 3) As Double 

‘thickness to cord ratio. 

Public k As Integer 

Dim o As Integer, p As Single 

Public Xvar() As Double 

Public Yvar() As Double 

Public Sub Incidence(i As Integer, j As Integer, slopegraph As Boolean, ksh As Single) 
‘Variation of (10)10-(delta0)10 
‘Load constants 
"solidity = 0.4 
i0delta0(6, 1) = -1.72244891587855E-11 
i0delta0(5, 1) = 3.05183134666209E-09 
i0delta0(4, 1) = -2.54582401992831E-07 
i0delta0(3, 1) = 8.9308266844057E-06 
i0delta0(2, 1) = -1.64583732868095E-04 
i0delta0(1, 1) = 2.36162560875073E-02 
" solidity = 0.6 
i0delta0(6, 2) = -1.21461996934098E- 10 
i0delta0(5, 2) = 2.44632274559037E-08 
i0delta0(4, 2) = -1.89083802948353E-06 
i0delta0(3, 2) = 6.51024967552871E-05 
i0delta0(2, 2) = -1.10036026944726E-03 
i0delta0(1, 2) = 4.37180120798075E-02 
‘ solidity = 0.8 
iddelta0(6, 3) = -1.56214579635869E-10 
i0delta0(5, 3) = 3.09609686430234E-08 
i0delta0(4, 3) = -2.3589555206982E-06 
i0delta0(3, 3) = 7.920583 15148836E-05 
i0delta0(2, 3) = -1.27482005365209E-03 
i0delta0(1, 3) = 5.69282884198401E-02 
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" solidity = 1.0 
i0delta0(6, 4) = -6.93736034940097E- I 1 
i0delta0(5, 4) = 8.79693724809005E-09 
i0delta0(4, 4) = -3.60272149557694E-07 
i0delta0(3, 4) = -4.42503626629787E-07 
i0delta0(2, 4) = 1.21654421718631E-05 
i0delta0(1, 4) = 6.52732436619772E-02 
" solidity = 1.2 
i0delta0(6, 5) = 1.04197027145803E- 10 
i0delta0(5, 5) = -2.10810745440021E-08 
i0delta0(4, 5) = 1.37089527130208E-06 
i0delta0(3, 5) = -3.80616685049517E-05 
i0delta0(2, 5) = 8.3487583 1969839E-05 
i0delta0(1, 5) = 8.43313899049463E-02 
‘ solidity = 1.4 
i0delta0(6, 6) = 5.1784789035559 1E-11 
i0delta0(5, 6) = -1.11055450426056E-08 
i0delta0(4, 6) = 6.3411319506379E-07 
i0delta0(3, 6) = -1.51366527276764E-05 
i0delta0(2, 6) = -1.39016582579643E-04 
i0delta0(1, 6) = 9.70746619623242E-02 
" solidity = 1.6 
i0delta0(6, 7) = -2.0832321334836 1E-10 
i0delta0(5, 7) = 3.637575297921 19E-08 
i0delta0(4, 7) = -2.51178491095239E-06 
i0delta0(3, 7) = 7.46391658452694E-05 
i0delta0(2, 7) = -1.21453189910881E-03 
i0delta0(1, 7) = 0.11686046293471 
" solidity = 1.8 
i0delta0(6, 8) = -1.56143738685847E-10 
iMdelta0(5, 8) = 2.10079058766965E-08 
i0delta0(4, 8) = -9.83116558606056E-07 
i0delta0(3, 8) = 6.88832778905635E-06 
i0delta0(2, 8) = 6.30398362773121E-05 
i0delta0(1, 8) = 0.124576532031824 
" solidity = 2 
i0delta0(6, 9) = -3.47201982252689E-10 
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i0delta0(5, 9) = 6.02515059199005E-08 
i0delta0(4, 9) = -4.22495870910922E-06 
i0delta0(3, 9) = 1.38861665249124E-04 
i0delta0(2, 9) = -2.50961937126704E-03 
i0delta0(1, 9) = 0.157387592875239 
k=9 

ReDim curve(k) 

ReDim Xvar(k) 

ReDim Yvar(k) 

p=04 

For o=1 Tok 


Yvar(o) = CurveFit2(i0delta0(1, 0), i0delta0(2, 0), iddelta0(3, 0), i0delta0(4, 0), 
10delta0(5, 0), 10delta0(6, 0), beta(i, j)) 


Xvar(o) =p 
p=p+0.2 

Next o 

QuadCoeff k 

10d0(i, j) = Interp(sigma(i, j)) 

If slopegraph = True Then 
"NACA 65-(A10)-series blades as equivalent circular arc 
‘Load constants 
" solidity = 0.4 
slopemn(7, 1) = 2.08333333236381E-11 
slopemn(6, 1) = -4.31089743405371E-09 
slopemn(5, 1) = 3.42948717823449E-07 
slopemn(4, 1) = -1.2816870624599E-05 
slopemn(3, 1) = 1.479268649298 1 5E-04 
slopemn(2, 1) = -3.97852565922108E-03 
slopemn(1, 1) = 0.53497377639842 
" solidity = 0.6 
slopemn(7, 2) = -3.47222221982991E-12 
slopemn(6, 2) = 7.772435880668E- 10 
slopemn(5, 2) = -5.87606836521815E-08 
slopemn(4, 2) = 1.3414189909966E-06 
slopemn(3, 2) = -5.4709838117617E-05 
slopemn(2, 2) = -2.55785261049368E-03 
slopemn(1, 2) = 0.67502768 1002291 
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" solidity = 0.8 

slopemn(7, 3) = -5.0821976835258E-21 
slopemn(6, 3) = -1.28205128154327E-10 
slopemn(5, 3) = 1.20192307773159E-08 
slopemn(4, 3) = -6.92016318915023E-07 
slopemn(3, 3) = -3.31075165860284E-05 
slopemn(2, 3) = -2.12211544896945E-03 
slopemn(1, 3) = 0.750014569226494 

" solidity = 1.0 

slopemn(7, 4) = -6.94444444304795E-12 
slopemn(6, 4) = 1.20192307712444E-09 
slopemn(5, 4) = -7.6655982672591E-08 
slopemn(4, 4) = 1.48528553811644E-06 
slopemn(3, 4) = -4.01670539815768E-05 
slopemn(2, 4) = -1.64423084743248E-03 
slopemn(1, 4) = 0.794953380497141 

" solidity = 1.2 

slopemn(7, 5) = -2.77777777637215E-11 
slopemn(6, 5) = 5.48076922735063E-09 
slopemn(5, 5) = -4.03579059504722E-07 
slopemn(4, 5) = 1.27039626889314E-05 
slopemn(3, 5) = -2.01934244898894E-04 
slopemn(2, 5) = -3.78525730411639E-04 
slopemn(1, 5) = 0.824973776832628 

' solidity = 1.4 

slopemn(7, 6) = 6.94444443457762E- 12 
slopemn(6, 6) = -9.1346153644617E- 10 
slopemn(5, 6) = 1.5758547106115E-08 
slopemn(4, 6) = 1.08027389256193E-06 
slopemn(3, 6) = -7.76942489579824E-05 
slopemn(2, 6) = -3.57051379069162E-04 
slopemn(1, 6) = 0.850023310682502 

" solidity = 1.6 

slopemn(7, 7) = -1.38888888844018E-11 
slopemn(6, 7) = 3.07692307570384E-09 
slopemn(5, 7) = -2.60683760439084E-07 
slopemn(4, 7) = 9.44930068769168E-06 
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slopemn(3, 7) = -1.87603922285007E-04 
slopemn(2, 7) = 6.39102459899732E-04 
slopemn(1, 7) = 0.870029 138228418 

" solidity = 1.8 

slopemn(7, 8) = -6.94444444643608E-12 
slopemn(6, 8) = 1.84294871811291E-09 
slopemn(5, 8) = -1.7841880317615E-07 
slopemn(4, 8) = 6.76354894046938E-06 
slopemn(3, 8) = -1.36371890420151E-04 
slopemn(2, 8) = 3.83012708653041E-04 
slopemn(1, 8) = 0.889956294448098 

" solidity = 2 

slopemn(7, 9) = -1.04166666679601E-11 
slopemn(6, 9) = 2.3637820511721E-09 
slopemn(5, 9) = -2.13141025329211E-07 
slopemn(4, 9) = 8.38760196941735E-06 
slopemn(3, 9) = -1.82296764080547E-04 
slopemn(2, 9) = 1.1642627027868E-03 
slopemn(1, 9) = 0.900013112661114 


k=9 
ReDim curve(k) 
ReDim Xvar(k) 
ReDim Yvar(k) 
p=04 
For o=1 Tok 
Yvar(o) = CurveFitl(slopemn(1, 0), slopemn(2, 0), slopemn(3, 0), 
slopemn(4, 0), slopemn(5, 0), slopemn(6, 0), slopemn(7, 0), 
beta(i, j)) 
Xvar(o) =p 
p=p+0.2 
Next o 
QuadCoeff k 


onemn(i, j) = Interp(sigmai, j)) 
Else 
‘Circular arc mean line blades 
‘Load constants 
" solidity = 0.4 
slopemn(7, 1) = -4.1666666643941 1E-11 
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slopemn(6, 1) = 8.14102563624798E-09 
slopemn(5, 1) = -5.809294868 14364E-07 
slopemn(4, 1) = 1.81949300572665E-05 
slopemn(3, 1) = -2.4787150325345E-04 
slopemn(2, 1) = -3.33942308314938E-03 
slopemn(1, 1) = 0.350002913822904 
" solidity = 0.6 
slopemn(7, 2) = 6.94444443965982E- 12 
slopemn(6, 2) = -1.4903846145501E-09 
slopemn(5, 2) = 1.06303418803688E-07 
slopemn(4, 2) = -3.08493589651349E-06 
slopemn(3, 2) = -5.89209358992093E-06 
slopemn(2, 2) = -2.75032054841517E-03 
slopemn(1, 2) = 0.54503205 1554482 
‘ solidity = 0.8 
slopemn(7, 3) = -1.04166666603368E-1 1 
slopemn(6, 3) = 1.56249999787653E-09 
slopemn(5, 3) = -7.29166664409364E-08 
slopemn(4, 3) = 1.18371202617595E-07 
slopemn(3, 3) = 2.38731068833431E-05 
slopemn(2, 3) = -3.17708338894818E-03 
slopemn(1, 3) = 0.649981061002052 
" solidity = 1.0 
slopemn(7, 4) = -6.94444443796575E-12 
slopemn(6, 4) = 7.21153844973621E-10 
slopemn(5, 4) = -2.93803398276893E-09 
slopemn(4, 4) = -2.18458625589335E-06 
slopemn(3, 4) = 3.707653560525 16E-05 
slopemn(2, 4) = -2.16987186519191E-03 
slopemn(1, 4) = 0.715017483006385 
' solidity = 1.2 
slopemn(7, 5) = 3.47222221644178E-12 
slopemn(6, 5) = -1.32211538505372E-09 
slopemn(5, 5) = 1.4369658127289E-07 
slopemn(4, 5) = -6.95767774239 187E-06 
slopemn(3, 5) = 1.06356110350703E-04 
slopemn(2, 5) = -2.01842956960263E-03 
7p 


slopemn(1, 5) = 0.760024767461402 

' solidity = 1.4 

slopemn(7, 6) = -6.7762635780344E- 21 
slopemn(6, 6) = -1.2820512967221E-10 
slopemn(5, 6) = 1.20192308605827E-08 
slopemn(4, 6) = -8.43531474004067E-07 
slopemn(3, 6) = -1.71984251977619E-05 
slopemn(2, 6) = -7.05448809185327E-04 
slopemn(1, 6) = 0.79493881180528 1 

‘ solidity = 1.6 

slopemn(7, 7) = 1.73611111042317E-11 
slopemn(6, 7) = -3.34134615158702E-09 
slopemn(5, 7) = 2.297008545038E-07 
slopemn(4, 7) = -7.7196241337063E-06 
slopemn(3, 7) = 9.70923186827122E-05 
slopemn(2, 7) = -1.42612189358715E-03 
slopemn(1, 7) = 0.825036422577128 
‘solidity = 1.8 

slopemn(7, 8) = 1.73611110991495E-11 
slopemn(6, 8) = -3.37339743378823E-09 
slopemn(5, 8) = 2.35309829010877E-07 
slopemn(4, 8) = -7.93524184317107E-06 
slopemn(3, 8) = 9.22482045950801E-05 
slopemn(2, 8) = -8.26442413767836E-04 
slopemn(1, 8) = 0.844992716324796 

' solidity = 2 

slopemn(7, 9) = 3.47222221135958E-12 
slopemn(6, 9) = -5.84935897461614E-10 
slopemn(5, 9) = 3.55235043048019E-08 
slopemn(4, 9) = -2.036349075 14072E-06 
slopemn(3, 9) = 3.84146771352789E-05 
slopemn(2, 9) = -6.86057805069 14E-04 
slopemn(1, 9) = 0.859969406337427 
k=9 

ReDim curve(k) 

ReDim Xvar(k) 

ReDim Yvar(k) 
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p=0.4 
For o= 1 Tok 


Yvar(o) = CurveFitl(slopemn(1, 0), slopemn(2, 0), slopemn(3, 0), 
slopemn(4, 0), slopemn(5, 0), slopemn(6, 0), slopemn(7, 0), 
beta(i, j)) 


Xvar(o) =p 
p=p+0.2 

Next o 

QuadCoeff k 

onemn(i, j) = Interp(sigmai, j)) 
End If 
' Variation of thickness-correction factor Kt for camber calculation 
‘Load constants 
‘betal = 10 
Kit(7, 1) = 816993.464355469 
Kit(6, 1) = -351150.076196289 
Kit(5, 1) = 61737.9966625977 
Kit(4, 1) = -5315.17809592133 
Kit(3, 1) = 142.41942977427 
Kit(2, 1) = 13.4581061812941 
Kit(1, 1) = 6.50000014506098E-02 
‘betal = 20 
Kit(7, 2) = 1077614.37939453 
Kit(6, 2) = -459489.065 126953 
Kit(5, 2) = 79133.1385789185 
Kit(4, 2) = -6710.73264049728 
Kit(3, 2) = 201.572413791658 
Kit(2, 2) = 12.2108848959033 
Kit(1, 2) = 7.74781832012336E-02 
betal = 30 
Kit(7, 3) = 796568.628417969 
Kit(6, 3) = -338565.234250488 
Kit(5, 3) = 58401.6780760498 
Kit(4, 3) = -4855.18434412186 
Kit(3, 3) = 103.6798390508 16 
Kit(2, 3) = 15.0572606819245 
Kit(1, 3) = 6.05909 105273496E-02 


‘betal = 40 
74 


Kit(7, 4) = -183823.529541016 
Kit(6, 4) = 137302.036435547 
Kit(5, 4) = -32831.3536986084 
Kit(4, 4) = 3879.52703175537 
Kit(3, 4) = -335.438802558397 
Kit(2, 4) = 26.0388210932938 
Kit(1, 4) = -3.53636350253408E-02 
‘betal = 50 
Kit(7, 5) = -449346.40869 1406 
Kit(6, 5) = 260840.875805664 
Kit(5, 5) = -56572.084102478 
Kit(4, 5) = 6297.7469984 1156 
Kit(3, 5) = -479.565571368521 
Kit(2, 5) = 30.7062532322958 
Kit(1, 5) = -7.44999987640929E-02 
‘betal = 60 
Kit(7, 6) = 2879901.96362305 
Kit(6, 6) = -1181513.95310303 
Kit(5, 6) = 193177.319282959 
Kit(4, 6) = -15656.7577867749 
Kit(3, 6) = 525.58354015793 
Kit(2, 6) = 8.91909924834078 
Kit(1, 6) = 0.1269090924257 12 
‘betal = 70 
Kit(7, 7) = 5187908.49780273 
Kit(6, 7) = -1994626.6972876 
Kit(5, 7) = 298546.694419434 
Kit(4, 7) = -21587.099794791 
Kit(3, 7) = 600.425619817768 
Kit(2, 7) = 13.8798799771508 
Kit(1, 7) = 9.83181832073486E-02 
k=7 
ReDim curve(k) 
ReDim Xvar(k) 
ReDim Yvar(k) 
p=10 
For o= 1 Tok 

iis) 


Yvar(o) = CurveFit1(Kit(1, 0), Kit(2, 0), Kit(3, 0), Kit(4, 0), Kit(5, 0), Kit(6, 0), 
Kit(7, 0), te(i, j)) 


Xvar(o) =p 

p=p+10 
Next o 
QuadCoeff k 
ktbar(i, j) = Interp(beta(i, j)) 
camber(i, j) = camberfunc(beta(i, j), beta((i - 1), j), ksh, ktbar(, j), i0d0G, j), onemn(i, j)) 
‘Zero-camber incidence angle 
‘Load constants 
" solidity = 0.4 
10(6, 1) = -4.04805786872846E- 14 
10(5, 1) = 3.30321577382553E-10 
10(4, 1) = -1.61193689951489E-07 
10(3, 1) = 1.60266684012811E-05 
10(2, 1) = -6.12705630658184E-04 
101, 1) = 3.96414794354314E-02 
" solidity = 0.6 
10(6, 2) = 2.77737297393952E- 10 
10(5, 2) = -5.8003011805044E-08 
10(4, 2) = 4.49158409132622E-06 
10(3, 2) = -1.58973331821244F-04 
10(2, 2) = 2.41923882140327E-03 
101, 2) = 3.77248 1264741 15E-02 
" solidity = 0.8 
10(6, 3) = -2.77737297421058E- 10 
10(5, 3) = 5.3836345 1353663E-08 
10(4, 3) = -3.86658408690899E-06 
10(3, 3) = 1.2355666447661E-04 
10(2, 3) = -1.73173878386024E-03 
i0(1, 3) = 7.33585 199368463E-02 
" solidity = 1.0 
10(6, 4) = -6.24048706766443E- 10 
10(5, 4) = 1.33904109689276E-07 
10(4, 4) = -1.07952815886492E-05 
10(3, 4) = 3.96289954 153417E-04 
10(2, 4) = -6.38683407 169083E-03 
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10(1, 4) = 0.11471689445807 
‘ solidity = 1.2 
10(6, 5) = 1.3907105163107E-10 
10(5, 5) = -2.856978047 1943 1E-08 
10(4, 5) = 2.01009383982154E-06 
10(3, 5) = -5.64950093888683E-05 
10(2, 5) = 3.77314309 162102E-04 
i0(1, 5) = 0.100321674263796 
‘ solidity = 1.4 
10(6, 6) = -1.04116066001981E-10 
10(5, 6) = 2.04204313875023E-08 
10(4, 6) = -1.67350788249365E-06 
10(3, 6) = 6.850833028693 16E-05 
10(2, 6) = -1.4205762442998 1E-03 
i0(1, 6) = 0.123885649683871 
" solidity = 1.6 
10(6, 7) = -6.59874024969446E- 10 
10(5, 7) = 1.26238706023296E-07 
10(4, 7) = -9.07669856742288E-06 
10(3, 7) = 3.00725006795233E-04 
10(2, 7) = -4.56882667 162972E-03 
i0(1, 7) = 0.153343047833914 
" solidity = 1.8 
10(6, 8) = -2.42843032737926E- 10 
10(5, 8) = 3.99324784211563E-08 
10(4, 8) = -2.4176220057015E-06 
10(3, 8) = 6.48183223432852E-05 
10(2, 8) = -8.29475889986497E-04 
i0(1, 8) = 0.149069730896372 
‘ solidity = 2 
10(6, 9) = -6.5928705639891E-10 
10(5, 9) = 1.27699043138418E-07 
10(4, 9) = -9.55189003892798E-06 
10(3, 9) = 3.42296645044371E-04 
10(2, 9) = -5.90334486059874E-03 
i0(1, 9) = 0.198583259410952 
k=9 
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ReDim curve(k) 
ReDim Xvar(k) 
ReDim Yvar(k) 
p=04 
For o=1 Tok 

Yvar(o) = CurveFit2(i0(1, 0), i0(2, 0), 10(3, 0), 10(4, 0), 10(5, 0), 10(6, 0), beta(i, 

)) 

Xvar(o) =p 

p=pt+0.2 
Ne xt o 
QuadCoeff k 
1010, j) = Interp(sigma(i, j)) 
‘Minimum loss incidence angle slope factor 
' Load constants 
" solidity = 0.4 
slopen(7, 1) = 5.20833333567409E- 12 
slopen(6, 1) = -1.00560897412991E-09 
slopen(5, 1) = 7.05128204714356E-08 
slopen(4, 1) = -2.35449445540326E-06 
slopen(3, 1) = 7.51638881979488E-06 
slopen(2, 1) = -3.64078518850874E-03 
slopen(1, 1) = -5.00123838208282E-02 
" solidity = 0.6 
slopen(7, 2) = -1.73611110737386E-12 
slopen(6, 2) = 2.604166663 12754E- 10 
slopen(5, 2) = -2.25694444522251E-08 
slopen(4, 2) = 1.0866477326843E-06 
slopen(3, 2) = -5.59643317217251E-05 
slopen(2, 2) = -2.28645827399987E-03 
slopen(1, 2) = -4.49905306238207E-02 
solidity = 0.8 
slopen(7, 3) = 6.94444444474202E- 12 
slopen(6, 3) = -1.85897435877984E-09 
slopen(5, 3) = 1.70806623789321E-07 
slopen(4, 3) = -7.15544871177087E-06 
slopen(3, 3) = 1.05446046163138E-04 
slopen(2, 3) = -2.7290063578107E-03 
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slopen(1, 3) = -4.00160259214886E-02 
" solidity = 1.0 

slopen(7, 4) = -5.2083333280508E-12 
slopen(6, 4) = 9.57532051044929E- 10 
slopen(5, 4) = -7.77243589489274E-08 
slopen(4, 4) = 3.08220425893069E-06 
slopen(3, 4) = -8.76518800509984E-05 
slopen(2, 4) = -8.65945466387075E-04 
slopen(1, 4) = -3.49963580541726E-02 
" solidity = 1.2 

slopen(7, 5) = 2.60416666588887E- 11 
slopen(6, 5) = -5.73317307430571E-09 
slopen(5, 5) = 4.60336538235517E-07 
slopen(4, 5) = - 1.68285620523179E-05 
slopen(3, 5) = 2.38823935717392E-04 
slopen(2, 5) = -2.12556085602955E-03 
slopen(1, 5) = -2.99817892565741E-02 
" solidity = 1.4 

slopen(7, 6) = 1.56249999943 168E-11 
slopen(6, 6) = -3.43349358927075E-09 
slopen(5, 6) = 2.74038461367532E-07 
slopen(4, 6) = -1.01127258078648E-05 
slopen(3, 6) = 1.33344623499454E-04 
slopen(2, 6) = - 1.06402239998715E-03 
slopen(1, 6) = -2.49992717449814E-02 
" solidity = 1.6 

slopen(7, 7) = 1.56249999951638E- 11 
slopen(6, 7) = -3.32131410091599E-09 
slopen(5, 7) = 2.54407051133998E-07 
slopen(4, 7) = -8.82776077215652E-06 
slopen(3, 7) = 9.46172051285998E-05 
slopen(2, 7) = -2.87900609066583E-04 
slopen(1, 7) = -1.99978148450413E-02 
" solidity = 1.8 

slopen(7, 8) = -5.20833332889783E-12 
slopen(6, 8) = 1.11778846096679E-09 
slopen(5, 8) = -9.01442307604805E-08 
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slopen(4, 8) = 2.80612616787579E-06 
slopen(3, 8) = -5.87438087080727E-05 
slopen(2, 8) = 3.75240412154199E-04 
slopen(1, 8) = -0.014986159811361 
"solidity = 2 

slopen(7, 9) = 1.21527777732163E-11 
slopen(6, 9) = -2.31971153754913E-09 
slopen(5, 9) = 1.6159188022391E-07 
slopen(4, 9) = -5.68345716178698E-06 
slopen(3, 9) = 7.84089690597511E-05 
slopen(2, 9) = -3.66426257727426E-04 
slopen(1, 9) = -9.98615979092321E-03 
k=9 

ReDim curve(k) 

ReDim Xvar(k) 

ReDim Yvar(k) 

p=04 

For o=1 Tok 


Yvar(o) = CurveFitl(slopen(1, 0), slopen(2, 0), slopen(3, 0), slopen(4, 0), 
slopen(5, 0), slopen(6, 0), slopen(7, 0), beta(i, j)) 


Xvar(o) =p 
p=p+0.2 
Next o 
QuadCoeff k 
NG, j) = Interp(sigma(i, j)) 
"Maximum thickness correction factor 
‘Load constants 
ikit(6) = -748795.36578369 1 
ikit(S) = 243951.67288208 
ikit(4) = -28087.2979736328 
ikit(3) = 1612.26135325431 
ikit(2) = -137.429116554558 
ikit(1) = 18.8187430176185 
ikt@, j) = CurveFit2(ikit(1), ikit(2), ikit(3), ikit(4), ikit(5S), ikit(6), tc(, j)) 
i0ref(i, j) = ksh * ikt(, j) * i010, j) 
i2d = i0ref(i, }) + NG, j) * camber(i, j) 
If i= 2 Then 
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If ksh = 0.7 Then 
icor(i, j) = 0.7238 * Mw(i- 1, j) + 7.5481 + 12d 
Elself ksh = 1.1 Then 
icor(i, j) = 1.3026 * Mw(i - 1, j) +. 5.738 + i2d 
Elself ksh = 1 Then 
icor(i, j) = i2d 
End If 
ElseIf i = 3 Then 
If ksh = 0.7 Then 
icor(i, j) = 0.7238 * MG - 1, j) + 7.5481 + 12d 
Elself ksh = 1.1 Then 
icor(i, j) = 1.3026 * MG - 1, j) + 5.738 + 12d 
Elself ksh = 1 Then 
icor(i, j) = i2d 
End If 
End If 
'Zero-camber deviation angle 
‘Load constants 
" solidity = 0.4 
d0(6, 1) = -6.95861265439764E- 11 
d0(5, 1) = 1.57394588647108E-08 
d0(4, 1) = -1.36278903151155E-06 
d0(3, 1) = 5.7135006738207E-05 
d0(2, 1) = -1.06391418466956E-03 
d0(1, 1) = 1.59951786608872E-02 
" solidity = 0.6 
d0(6, 2) = -7.08410313374729E- 14 
d0(5, 2) = -4.63603904397869E- 10 
d0(4, 2) = 1.86661040402214E-07 
d0(3, 2) = -1.62241633923088E-05 
d0(2, 2) = 6.30890134971196E-04 
d0(1, 2) = 5.06008922957335E-03 
" solidity = 0.8 
d0(6, 3) = -1.72852100539214E- 11 
d0(5, 3) = 3.54731371186856E-09 
d0(4, 3) = -1.83872938541718E-07 
d0(3, 3) = 1.7208295339799E-06 
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d0(2, 3) = 2.60107809708643E-04 
d0(1, 3) = 9.95347546358971E-03 
" solidity = 1.0 

d0(6, 4) = 1.91189805496704E- 10 
d0(5, 4) = -3.92754785100841E-08 
d0(4, 4) = 3.08863831766093E-06 
d0(3, 4) = -1.10622510739233E-04 
d0(2, 4) = 2.00666089 165225E-03 
d0(1, 4) = 9.999633 12096952E-04 
' solidity = 1.2 

d0(6, 5) = -2.95336231991993E- 10 
d0(5, 5) = 6.09853 177350322E-08 
d0(4, 5) = -4.62262480382947E-06 
d0(3, 5) = 1.62505009349445E-04 
d0(2, 5) = -2.2878080447981 1E-03 
d0(1, 5) = 2.54709629980425E-02 
' solidity = 1.4 

d0(6, 6) = -1.75483338526582E- 11 
d0(5, 6) = 1.5277372943 101 1E-09 
d0(4, 6) = 2.26701400096729E-07 
d0(3, 6) = -2.11891579624535E-05 
d0(2, 6) = 8.81687814285215E-04 
d0(1, 6) = 7.87309303404982E-03 
" solidity = 1.6 

d0(6, 7) = -3.46716215172916E- 11 
d0(5, 7) = 7.92043 135938725E-09 
d0(4, 7) = -5.62396773595708E-07 
d0(3, 7) = 2.2674997246952E-05 
d0(2, 7) = -1.57381818780777E-04 
d0(1, 7) = 0.017343983408864 

‘ solidity = 1.8 

d0(6, 8) = 6.30485119950015E-11 
d0(5, 8) = -1.41425240546278E-08 
d0(4, 8) = 1.37167452307629E-06 
d0(3, 8) = -5.62863600421792E-05 
d0(2, 8) = 1.32778646161569E-03 
d0(1, 8) = 8.46711704980407E-03 
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"solidity = 2 

d0(6, 9) = 8.67296545047297E- 1 1 
d0(5, 9) = -1.60473137209016E-08 
d0(4, 9) = 1.29498405332384E-06 
d0(3, 9) = -4.33874967669112E-05 
d0(2, 9) = 8.78086666489253E-04 
d0(1, 9) = 1.43381905636488E-02 
k=9 

ReDim curve(k) 

ReDim Xvar(k) 

ReDim Yvar(k) 

p=04 

Foro=1 Tok 


Yvar(o) = CurveFit2(d0(1, 0), d0(2, 0), d0(3, 0), d0(4, 0), d0(5, 0), d0(6, 0), 
beta(i, j)) 


Xvar(o) =p 
p=pt+0.2 

Next o 

QuadCoeff k 

d010G, j) = Interp(sigma(i, j)) 

If slopegraph = True Then 
"NACA 65-(A10)-series blades as equivalent circular arc 
‘Load constants 
" solidity = 0.4 
slopem(7, 1) = 6.94444443457762E-12 
slopem(6, 1) = -1.49038461433326E-09 
slopem(5, 1) = 1.06303418845322E-07 
slopem(4, 1) = -3.16069348027526E-06 
slopem(3, 1) = 7.70624526893471E-05 
slopem(2, 1) = -7.91987265529315E-04 
slopem(1, 1) = 0.412494173012831 
" solidity = 0.6 
slopem(7, 2) = 2.08333333206735E-1 1 
slopem(6, 2) = -4.48717948388355E-09 
slopem(5, 2) = 3.5296474354407 1E-07 
slopem(4, 2) = -1.22224650365155E-05 
slopem(3, 2) = 2.05185752662373E-04 
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slopem(2, 2) = -8.88621860667627E-04 
slopem(1, 2) = 0.2774985435 14405 

' solidity = 0.8 

slopem(7, 3) = 1.04166666611838E-11 
slopem(6, 3) = -2.20352563908 184E-09 
slopem(5, 3) = 1.64262820512295E-07 
slopem(4, 3) = -4.77163461454211E-06 
slopem(3, 3) = 6.08733981977139E-05 
slopem( 2, 3) = 2.12339688573593E-04 
slopem(1, 3) = 0.210016025961522 

" solidity = 1.0 

slopem(7, 4) = 3.47222221813584E-12 
slopem(6, 4) = -1.19391025581519E-09 
slopem(5, 4) = 1.31677350481696E-07 
slopem(4, 4) = -5.84899476141487E-06 
slopem(3, 4) = 1.26963627380405E-04 
slopem(2, 4) = -5.62980817107928E-04 
slopem(1, 4) = 0.170010198410182 

" solidity = 1.2 

slopem(7, 5) = 5.2083333305919E-12 
slopem(6, 5) = -1.51842948674721E-09 
slopem(5, 5) = 1.550480769885 18E-07 
slopem(4, 5) = -6.59036276573488E-06 
slopem(3, 5) = 1.3337230542021E-04 
slopem(2, 5) = -5.56330170923047E-04 
slopem(1, 5) = 0.142489073668344 

‘ solidity = 1.4 

slopem(7, 6) = 2.25694444314355E-11 
slopem(6, 6) = -4.79567307436549E-09 
slopem(5, 6) = 3.83947649437721E-07 
slopem(4, 6) = -1.38976908479194E-05 
slopem( 3, 6) = 2.37974820819886E-04 
slopem(2, 6) = -1.0776442692304E-03 
slopem(1, 6) = 0.122499271780747 

" solidity = 1.6 

slopem(7, 7) =-1.04166666620308E- 11 
slopem(6, 7) = 1.94711538320687E-09 
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slopem(5, 7) = -1.29807692192396E-07 
slopem(4, 7) = 4.01260197557818E-06 

slopem(3, 7) = -4.06300984678865E-05 
slopem(2, 7) = 3.55929451188786E-04 

slopem(1, 7) = 0.107513112087279 

" solidity = 1.8 

slopem(7, 8) = 1.21527777736398E-11 

slopem(6, 8) = -2.80048076742312E-09 
slopem(5, 8) = 2.45726495706755E-07 

slopem(4, 8) = -9.71317744902223E-06 
slopem(3, 8) = 1.8301998394321 1E-04 

slopem(2, 8) = -8.29567340645099E-04 
slopem(1, 8) = 9.50021854990979E-02 

" solidity = 2 

slopem(7, 9) = -3.47222222152398E- 12 
slopem(6, 9) = 6.00961538020119E-10 

slopem(5, 9) = -3.31196580577453E-08 
slopem(4, 9) = 8.84688224900287E-07 

slopem(3, 9) = 4.18366859378239E-07 

slopem(2, 9) = 1.88301250858558E-04 

slopem(1, 9) = 8.74956295422606E-02 


k=9 
ReDim curve(k) 
ReDim Xvar(k) 
ReDim Yvar(k) 
p=04 
For o= 1 Tok 
Yvar(o) = CurveFitl(slopem(1, 0), slopem(2, 0), slopem(3, 0), 
slopem(4, 0), slopem(5, 0), slopem(6, 0), slopem(7, 0), beta(i, 
)) 
Xvar(o) =p 
p=p+0.2 
Next o 
QuadCoeff k 


dm(i, j) = Interp(sigma(i, j)) 
Else 
' Circular arc mean line blades 


‘Load constants 


85 


" solidity = 0.4 

slopem(7, 1) = -5.0821976835258E-21 
slopem(6, 1) = 6.73076922756005E-10 
slopem(5, 1) =-1.1258012813653E-07 
slopem(4, 1) = 6.10941141054866E-06 
slopem(3, 1) = -1.28908069711997E-04 
slopem(2, 1) = 1.04214733619301E-03 
slopem(1, 1) = 0.602489802490481 

' solidity = 0.6 

slopem(7, 2) = 5.20833332212157E-12 
slopem(6, 2) = -1.08573717811505E-09 
slopem(5, 2) = 8.45352564893265E-08 
slopem(4, 2) = -3.08293269490889E-06 
slopem(3, 2) = 5.90424689903557E-05 
slopem(2, 2) = -9.9919948297611E-05 
slopem(1, 2) = 0.409991987649079 

‘ solidity = 0.8 

slopem(7, 3) = 1.73611110229166E-12 
slopem(6, 3) = -1.9631410169349E- 10 
slopem(5, 3) = 9.34829047505303E-10 
slopem(4, 3) = 6.51405884788403E-07 
slopem(3, 3) = -1.73966821961358E-05 
slopem(2, 3) = 6.07932628 156505E-04 
slopem(1, 3) = 0.310002185699076 
"solidity = 1.0 

slopem(7, 4) = -6.94444444474202E- 12 
slopem(6, 4) = 1.66666666600625E-09 
slopem(5, 4) = -1.47569444275408E-07 
slopem(4, 4) = 6.10795453948043E-06 
slopem(3, 4) = -1.06196337540609E-04 
slopem(2, 4) = 1.0979166113998E-03 
slopem(1, 4) = 0.24998 1060934445 

" solidity = 1.2 

slopem(7, 5) = -1.21527777719458E- 11 
slopem(6, 5) = 2.57612179309884E-09 
slopem(5, 5) = -2.01255341752971E-07 
slopem(4, 5) = 7.285292825898 1 8E-06 
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slopem(3, 5) = -1.10063251923975E-04 
slopem(2, 5) = 1.03774033499349E-03 
slopem(1, 5) = 0.210013840616803 
" solidity = 1.4 
slopem(7, 6) = -6.94444444135388E- 12 
slopem(6, 6) = 1.53846153871928E-09 
slopem(5, 6) = -1.19925213509786E-07 
slopem(4, 6) = 4.09965034364745E-06 
slopem(3, 6) = -4.4843627868829 1E-05 
slopem(2, 6) = 3.82051237181713E-04 
slopem(1, 6) = 0.182514569024661 
" solidity = 1.6 
slopem(7, 7) = -2.25694444348236E- 11 
slopem(6, 7) = 4.50721153650614E-09 
slopem(5, 7) = -3.28258546761218E-07 
slopem(4, 7) = 1.07734192780207E-05 
slopem(3, 7) = -1.40615408163569E-04 
slopem(2, 7) = 8.60176240280452E-04 
slopem(1, 7) = 0.160005099306115 
" solidity = 1.8 
slopem(7, 8) = 1.73611110991495E-12 
slopem(6, 8) = -5.48878204605779E- 10 
slopem(5, 8) = 6.78418804189729E-08 
slopem(4, 8) = -3.58591929172647E-06 
slopem(3, 8) = 9.18654094732574E-05 
slopem(2, 8) = -3.68509653696947E-04 
slopem(1, 8) = 0.142494901151593 
‘ solidity = 2 
slopem(7, 9) = -5.20833333567409E- 12 
slopem(6, 9) = 1.14983974338484E-09 
slopem(5, 9) = -9.05448717253288E-08 
slopem(4, 9) = 3.42894084548462E-06 
slopem(3, 9) = -5.49887086620515E-05 
slopem(2, 9) = 7.85977528323656E-04 
slopem(1, 9) = 0.1275007286422 
k=9 
ReDim curve(k) 
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ReDim Xvar(k) 
ReDim Yvar(k) 
p=04 
For o= 1 Tok 
Yvar(o) = CurveFitl(slopem(1, 0), slopem(2, 0), slopem(3, 0), 
slopem(4, 0), slopem(5, 0), slopem(6, 0), slopem(7, 0), beta(i, 
)) 
Xvar(o) =p 
p=p+0.2 
Next o 
QuadCoeff k 
dm(i, j) = Interp(sigma(i, j)) 
End If 
"Maximum thickness correction factor 
‘Load constants 
dkit(6) = 618823.625244141 
dkit(5) = -202775.302703857 
dkit(4) = 25013.8597869873 
dkit(3) = -1269.01561832427 
dkit(2) = 41.3428950682282 
dkit(1) = 7.56794627627824 
dkt(i, j) = CurveFit2(ikit(1), ikit(2), ikit(3), ikit(4), ikit(S), ikit(6), tc(i, j)) 
dOref(i, j) = ksh * dkt(i, j) * d010Q, j) 
dref(i, j}) = dOref(i, j) + dmQ, j) * camber(i, j) 
End Sub 


14. Module mdlInterpolation Code 


Option Explicit 
Option Base | 
Public Ainterp(Q) As Double 
Public Binterp() As Double 
Public Cinterp() As Double 
Dim o As Integer 
Public Sub QuadCoeff(k As Integer) 
ReDim Ainterp(k) 
ReDim Binterp(k) 
ReDim Cinterp(k) 
Foro=2Tok- 1 
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Ainterp(o) = interpA(Xvar(o + 1), Xvar(o - 1), Xvar(o), Yvar(o), Yvar(o + 1), 
Yvar(o - 1)) 


Binterp(o0) = interpB(Yvar(o), Yvar(o - 1), Xvar(o), Xvar(o - 1), Ainterp(o)) 
Cinterp(o0) = interpC(Yvar(o), Ainterp(o), Binterp(o), Xvar(o)) 
Next o 
Ainterp(1) = Ainterp(2) 
Binterp(1) = Binterp(2) 
Cinterp(1) = Cinterp(2) 
Ainterp(k) = Ainterp(k - 1) 
Binterp(k) = Binterp(k - 1) 
Cinterp(k) = Cinterp(k - 1) 
Foro=1Tok- 1 
Ainterp(o) = (Ainterp(o) + Ainterp(o + 1))/2 
Binterp(o) = (Binterp(o) + Binterp(o + 1))/2 
Cinterp(o0) = (Cinterp(o) + Cinterp(o + 1)) /2 
Next o 
End Sub 
Public Function Interp(Xval As Double) 
o=2 
Do 
If Xval < Xvar(o) Then 
Interp = Ainterp(o) * Xval * 2 + Binterp(o) * Xval + Cinterp(o) 


Exit Do 
Else 
If Xval > Xvar(o) And Xval < Xvar(o + 1) Then 
Interp = Ainterp(o) * Xval * 2 + Binterp(o) * Xval + 
Cinterp(o) 
Exit Do 
Else 
o=o0+l 
End If 
End If 
Ifo=k-1 Then 
Interp = Ainterp(o) * Xval * 2 + Binterp(o) * Xval + Cinterp(o) 
Exit Do 
End If 


Loop Until o =k 


End Function 


89 


15. Module mdlFunctions Code 


Option Exp licit 
Option Base | 


Public Function Diffusion(phil As Double, betal As Double, beta2 As Double, pitch_ml As 
Double, solidity As Double) As Double 


Diffusion = 1 - phil * (Cos(DegToRad(betal)) / Cos(DegToRad(beta2))) + 
(((Tan(DegToRad(betal)) - pitch_ml * phil * Tan(DegToRad(beta2))) * 
Cos(DegToRad(betal))) / ((1 + pitch_ml) * solidity)) 


End Function 
Public Function ggl(g As Double) As Double 
ggl=g/(g- 1) 
End Function 
Public Function g1(g As Double) As Double 
gl=1/(g-1) 
End Function 
Public Function g1g(g As Double) As Double 
glp=(g- 1712 
End Function 
Public Function Mach(g As Double, X1 As Double) As Double 
Mach = Sqr(((2 / (g- 1)) * (K1% 2))/(1 - (KI % 2))) 
End Function 
Public Function DegToRad(angle As Double) As Double 
DegToRad = angle * ((22 /7) / 180) 
End Function 
Public Function RadToDeg(radians As Double) As Double 
RadToDeg = radians * (180 / (22 /7)) 
End Function 
Public Function Arccos(ratio As Double) As Double 
Arccos = RadToDeg(Atn(-ratio / Sqr(-ratio * ratio + 1)) + 2 * Atn(1)) 
End Function 
Public Function Arcsin(ratio As Double) As Double 
Arcsin = RadToDeg(Atn(ratio / Sqr(-ratio * ratio + 1))) 
End Function 


Public Function DofReaction(Xthetal As Double, Xul As Double, Xtheta2 As Double, XU2 As 
Double) As Double 


DofReaction = (1 - 0.5 * ((Xthetal / Xul) + (Xtheta2 / XU2))) 


End Function 
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Public Function SFTC(betal As Double, beta2 As Double, solidity As Double, height As Double, 
tipgap As Double, spacing As Double) 


Dim CDi As Double 

Dim CL As Double 

Dim betainf As Double 

betainf = Atn((Tan(DegToRad(betal)) + Tan(DegToRad(beta2))) / 2) 


CL = (2 / solidity) * (Tan(DegToRad(betal)) - Tan(DegToRad(beta2))) * 
Cos(DegToRad(betainf)) 


CDi = (0.25 * (CL * 2) * solidity * (tipgap / height) * (1 / Cos(DegToRad(beta2)))) + 
(0.04 * (CL “ 2) * solidity * (spacing / height)) 


SFTC = CDi * ((Cos(DegToRad(beta1)) * 2) / (Cos(DegToRad(betainf)) * 3)) * solidity 
End Function 


Public Function taufunc(XU2 As Double, Xtheta2 As Double, Xul As Double, Xthetal As 
Double) As Double 


taufunc = 1 + 2* (XU2 * Xtheta2) - (Xul * Xthetal)) 

End Function 

Public Function Xfunc(Xz1 As Double, alphal As Double) As Double 
Xfunc = Xz1 / Cos(DegToRad(alpha1)) 

End Function 

Public Function Xthetafunc(rl As Double, r2 As Double, Xthetal As Double) As Double 
Xthetafunc = (rl / r2) * Xthetal 

End Function 

Public Function Xufunc(rl As Double, r2 As Double, Xul As Double) As Double 
Xufunc = (r2 /rl) * Xul 

End Function 

Public Function Xwfunc(Xz1 As Double, betal As Double) As Double 
Xwfunc = Xz1 / Cos(DegToRad(beta1)) 

End Function 

Public Function Yfunc(X1 As Double, taul As Double) As Double 
Yfunc = X1 / Sqr(taul) 

End Function 

Public Function Machz(M1 As Double, alphal As Double) As Double 
Machz = M1 * Cos(DegToRad(alphal )) 

End Function 

Public Function alphafunc(Xthetal As Double, Xz1 As Double) As Double 
alphafunc = RadToDeg(Atn(Xthetal / Xz1)) 

End Function 


Public Function betafunc(Xul As Double, Xthetal As Double, Xz1 As Double) As Double 


91 


betafunc = RadToDeg(Atn((Xul - Xthetal) / Xz1)) 

End Function 

Public Function rhfuncl(rt As Double, Al As Double) As Double 
rhfuncl = Sqr(rt * 2 - (Al /(22/7))) 

End Function 

Public Function rhfunc2(rm As Double, rt As Double) As Double 
rhfunc2 = 2 * rm - rt 

End Function 

Public Function rtfunc(Al As Double, rm As Double) As Double 
rtfunc = (Al /(4 * rm * (22/7)))+rm 

End Function 

Public Function rmfunc(rt As Double, rh As Double) As Double 
rmfunc = (rt + rh) /2 

End Function 

Public Function rhtfunc(rt As Double, rh As Double) As Double 
rhtfunc = rh /rt 

End Function 

Public Function rhtfunc2(h As Double, rm As Double) As Double 
Dim hrm As Double 
hrm = h / (2 * rm) 
rhtfunc2 = (1 - hrm) / (1 + hrm) 

End Function 


Public Function CurveFitl(constO As Double, const! As Double, const2 As Double, const3 As 
Double, const4 As Double, const5 As Double, const6 As Double, polyvar As Double) 


"No y intercept 


CurveFitl = (const6 * (polyvar * 6)) + (const5 * (polyvar “ 5)) + (const4 * (polyvar * 4)) 
+ (const3 * (polyvar * 3)) + (const2 * (polyvar * 2)) + (constl * polyvar) + 
const0 


End Function 


Public Function CurveFit2(constl As Double, const2 As Double, const3 As Double, const4 As 
Double, const5 As Double, const6 As Double, polyvar As Double) 


‘y intercept 


CurveFit2 = (const6 * (polyvar * 6)) + (const5 * (polyvar “ 5)) + (const4 * (polyvar * 4)) 
+ (const3 * (polyvar * 3)) + (const2 * (polyvar * 2)) + (const! * polyvar) 


End Function 


Public Function camberfunc(beta2 As Double, betal As Double, kshape As Single, kthick As 
Double, idelta As Double, slope As Double) 


camberfunc = ((beta2 - betal) - (kshape * kthick * idelta)) / slope 


End Function 
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Public Function interpA(Xip1 As Double, Xim1 As Double, Xi As Double, Yi As Double, Yip1 
As Double, Yim1 As Double) 


interpA = (1 / (Xip1 - Xim1)) * (((Yip1 - Yi) / (Xip1 - Xi)) - ((Yi- Yim1) / (Xi - Xim1))) 
End Function 


Public Function interpB(Yi As Double, Yim1 As Double, Xi As Double, Xim1 As Double, Ai As 
Double) 


interpB = ((Yi- Yim1) / (Xi - Xim1)) - Ai * (Xi + Xim1) 

End Function 

Public Function interpC(Yi As Double, Ai As Double, Bi As Double, Xi As Double) 
interpC = Yi- Ai * Xi* 2- Bi* Xi 

End Function 

Public Function ShockLoss(M As Double, g As Double) 
Dim ystar As Double 


ystar = (1/(4* g*M%2))* ((g+1)*M%2)-(- g)+Sar((g +1) * (e+ *M%4) 
“(2% (37 8) M2) 4g 2))) 


ShockLoss = ((g + 1)/(2* g *M% 2 * ystar - (g - 1)) * gl(g) *((g+1)*M%2* 
ystar) / (2 + (g- 1) *M% 2 * ystar)) * ggl(g) 


End Function 

Public Function DiffB(sigma As Double, R As Double) 
DiffB = 2 * sigma * R * (1 + R) 

End Function 


Public Function DiffA(sigma As Double, R As Double, D As Double, phi As Double, betal As 
Double) 


DiffA = (1 + R) * sigma * (1 - D) / (phi * Cos(DegToRad(betal)))) + 
(Tan(DegToRad(beta1)) / phi) 


End Function 
Public Function SinB2(A As Double, B As Double, R As Double) 


SinB2 = (-B + Sqr(B *2+(4* (A%2-(B*2/(4*R%2))) *(RA2+A%2)))/(2*(R 
A2+A42)) 


End Function 
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APPENDIX C. CALCULATION RESULTS 
HAND CALCULATIONS 


Inputs derived from Sanger [Ref. 13]. 
m = 17.091bm/sec 


P, =14.69psia 
T, =517.8R 
M,,,, = 1.28 


= 2836.33rad/sec 
R =53.35ft- Ibsf /Ibm-’ R 
y=1.4 


U, =0-r, = 2836.33 - (5.5 inches/12 inches per ft) = 1300 ft/sec 





U y 1.4 
Xj, =F —"—— where C, =| — |. R=| —* }-53.35= 186.725 
eS CET, ' | a4) 
1 
6 a 299210 
/2-186.725-32.2-517.8 
y-1 . 2 14-1... 
sin’ B -( Xin ae Blas. -( 0.5210" wie Palaans 
no (ll apes] laine Aq engne ll =—faep oo... 
14+Xi,, My, 1+0.5210 = | ag? 


sin’ B,, = 0.8649 
sinB,, =0.9300 = B, = sin~'(0.9300) = 68.438 


M.,,, =M ,,,, -cosB,, =1.28-cos(68.438") = 0.4704 


Zit Wi 


D,,, = 0.455; R,, = 1.0687; 6,, = 0.956 


2 


G,,, = 1.52; f, (6, /6,,) = 0.8553 
8, = 0.0045; AR, =1.2 


(t/c) = 0.037 at tip; 0.056 at mean; 0.08 at hub 


2max 


D,,, = 0.52; R,, = 1.0259; ,, = 1.0371 


3 
0;, =1.25;f, (63 /6,,,) =0.8 
5, = 0.003; AR, =1.2 


(t/c) = 0.07 at tip; 0.06 at mean; 0.05 at hub 


3max 


Applying equations from Appendix A yields the followings results. 
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1. Inlet Conditions 


M = 0.4704 





= 0.4704 
cosO 





X,,, =(0.2059)-cos0 = 0.2059 
X,,, = (0.2059). tan0 =0 
X,,, = 0+(0.2509) -tan(68.438) = 0.5210 


_ 14.69 
"  53.35-517.8 


V,, =N2-32.2 -186.725-517.8 = 2495.3139 ft/sec 


-144 = 0.07658 Ibm/ft’ 


1 
®,, = 0.2059-(1- 0.2059" )F #1 = 0.1847 


17.09 
A, =| ——————_. |. —————_ = 0.484 ft’ or 69.7144 in’ 
0.07658 -2495.3139 } 0.1847-cos0 


_ 0.5210 -2495.3139 


te 048 tt oso 
2836.33 


1, =4/5.5° —09.714 “= 2.8389 in or 0.2366 ft 


1, = 28389 i 5 = 0.5162 





542. 
lim = ae 4.1695 in or 0.3475 ft 
@lm = = ; 0 = 0 
4.1695 


X,,, =X, = 0.2058 


35 
o.,,, = tan’'| ———- tan0 |=0 
4.1695 


sg = tO? .0,5210= 0.3949 
BS hS 


3949 — 
B,,, = tan! om) = 62.4685 





0.2058 


2 0:2058 _ 9 2958 


cosO 
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a 0.2058" 





Maite eo 4704 
1—0.2058" 
— 9.2058 _ 9 4454 


wim C0862.4685 





2. Rotor Conditions at Mean Line 


cos62.4685 es 
cos Bois 
tan62.4685 — 1.0687 -0.956 - tan B,,, 
(1+ 1.0687) -1.52 


0.455 =1-(0.956)- 
-c0s62.4685 


Solve for B 


2m 


—-B+ *+4(R, +A) AB 2. | 
. = 21 
sinB,,, = a(R? +A?) 


Where A = (1 +R,, )-6,,, || + tanB,,, 
>,, -cosB,,, o,, 


1-0.455 }: tan62.4685 


A =| (1+1.0687)-1.52-| —————>__ 
0.956 -cos62.4685 0.956 


} 5.8848 


and B =2-6,,,-R,,-(1+R,, )= 2-1.52 -1.0687 -(1+ 1.0687) = 6.7209 


2m 


6.7209° 
4 -1.0687° 





—6.7209+ ,{6.7209° +4(1.0687° +5.8848° ) | sssas" - 


sinB,,, = 
i. 2 -(1.0687° + 5.8848’) 


sinB,,, = 0.7430 
B,,, =47.9895° 


X om = 1.0687 -0.3949 =0.4221 


U2) 


_ 0.2059 
™ 0.3949 





= 0.5213 


1 
05, = 0.956 -0.5213-| ——— |= 0.4663 
1.0687 


X 4, = 0.4663 - 0.4213 = 0.1968 


Z2 


X yom = 0.4221 — 0.1968tan47.9895 = 0.2036 


62 
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0.1968 


wo[ed ( 0.0 vase Grane 


0.2036 
o,,. = tan” ( } 45.9677" 





20.3949 0.4221 
t= 1+2-[0.4221 0.2036— 0.3949-0.0] =1.1718 


0.1968 


om a = 02.2832 
cos(45.9677) 


0.1968 


wom = 0.2941 
cos(47.9894) 


as OMENS =H 6tb 


Y 
™  ALATIS 


_ 0.2941 


Vo eee 
wn 1.1718 


= 0.2717 





45.9677 
ice ps 

cos0 
(tan45.9677—1.0259- 1.0371-tan0) -cos45.9677 
(1+ 1.0259) -1.25 


D 
= 0.5631 


X,,,, =1.0371 0.1968 =0.2041 


Z3m 


X,,, =0.2041- tan0 =0 


63m 


X ,,, = 0.4221-1.0259 = 0.4330 


U3 


B,,, = tan” OTe | 647615" 
0.2041 


_ 0.2041 
cosO 


= 0.2041 





3m 


= 020M 20.4787 


cos64.7615 


W3m 
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0.2041 


Y,, == = 0.1886 
eof 1718 
JOA: seas 


Y 
we f.1718 





4. Iteration #1 
Rotor Performance 
Initially set @, and @.,,.. =0 


2 
 =2-1.52- cos’ 62.4685 


i : ; : [ 0.005 +0.16 0.455° | = 0.02569 
7 cos’ 47.9894 


O, = 0+ 0+ 0.02569 


yA =1- 0.2059" =0.9576 
tl 


P 3. 
vA = (0.9576) =0.8594 


tl 


Th = 0.9576 + 0.4454’ = 1.1560 


1 
Pu =(1.1560)° =1.6608 
6 = (1 =f 


Tf =1.1560+ 0.4221° —0.3949* =1.1781 
tl 


A = 1781 


Pp 3:5 
"% =(1.1781) =1.7750 


tl 


Pf, = 1.7750 — 0.02569 - [1.6608 — 0.8594] = 1.7544 


tl 


fe 
1.1718 
Puf =1.7544- =1.7218 
LA 1.1781 


T, T= 1.1718—0.2832* =1.0917 
tl 


3.5 
1.0917 
“p =1.7218.( 7) 1.3435 
j L.1718 
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2.5 
®,,, = 0.2616-(1-0.2616°)  =0.2191 
Ay ___0.1847-cos0 TITS <3 oo 
A, 0.2191-cos45.9677 1.7218 
A, = 0.7627 - 69.7144 = 53.1725 in? or 0.3693 ft? 
1, =4.1694- 1.0687 = 4.4559 in or 0.3713 ft 


(= 82S _ — 1 9092 in or 0.1583 ft 
2-n-4,4559 
1.8992 


___2-4.4559 _ 9 6487 


Tyo 1.8992 
1+————_ 
2-4,4559 


r,, =| ——— |-4.4559 = 5.4055 in or 0.4505 ft 
1+ 0.6487 


r,, =0.6487 -5.4055 = 3.5063 in or 0.2922 ft 


Stator Performance 
Initially set @, and Oya. = 0 


cos’ 45.9677 


@, =21.25. - -[0.005+ 0.16- 0.5631* | = 0.02546 


©, =0+0+0.02546 


‘4 = 1.7218 —0.02546- (1.7218-1.3435) =1.7121 


1 
Sf, =T=1.1718 


T, T =1.1718-—0.2041° =1.1302 


u 


1.1302)” 
¥ = = 1.5084 
. | 1.1718 


@,, = 0.1886-(1-0.1886?) » = 0.1722 





A _ 0.1847-cos0 V1.1718 ~ 0.6782 
1 


~ 0.1722-cosO 1.7121 
A, = 0.6782: 69.7144 = 47.2778 in? or 0.3233 ft? 


r,,, =1.0259 -4.4559 = 4.5713 in or 0.3809 ft 


3 


47.2778 


3 == = 1.6460 in or 0.1372 ft 
2-%-4.5713 


100 


1.6460 


Belo CL ey 

ths Leeoe 
jee 
2-4.5713 


2 
r,, =| ————— }|-4.5713 =5.3943 in or 0.4495 ft 
1+0.6949 
r,, = 0.6949. 5.3943 = 3.7483 in or 0.3124 ft 


Rotor Hub Calculations 


% ,./=02036 es 
3.5063 


=0.2587 





X,., = 0.4221 22083 9.3391 
4.4559 


2587 
a, =ta0"(5 as } 52.7381" 





0.1968 
3321—0.2587 
B,, = tan”! O8e2T D287 | 504560" 
0.1968 
= 0.1968 anys 
cos52.7381 
or ~_ 0.1968 _ 9 540) 
cos20.4569 
0.3251 
Y,, =. = 0.3003 
eS ATIS. 
0.2101 
Yu = === = 0.1940 
m™ LA718 





Rotor Tip Calculations 


X,,, = 0.2036 4°? 0.1678 
5.4055 


Ko 0d PO? 
4.4559 





= 0.5120 
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0.1678 
o,, = tan” = 40.4526" 
0.1968 


ee 0.5120—0.1678 
a 0.1968 


- 60.2402’ 


0.1968 


oy = —@—— = 0.2586 
co0s40.4526 


_ 0.1968 _ 
Y  60860.2402 


3965 


Y= 9.2586 _ 0.2389 


* AL.1718 


w2t 








My 
ae 45713 _ 99 
3.7483 
Ho), 043302 94580 
; 5713 





a aa 0.3550—0.0 
ss 0.2041 


} 60.1058" 


_ 0.2041 
cosO 


Xx 


3h 





= 0.2041 


— 0.2041 
“" ¢0860.1058 


ya OM 5 1886 


“1.1718 


0.4095 
Yuya, === = 0.3783 


1.1718 


= 0.4095 
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_* 9.18867 
14-1 


1—0.1886° 






= 0.4293 


ea 0.3783" 
14-1 


M = 4, 
a 1— 0.37837 


= 0.9139 


Stator Tip Calculations 


x, =0.0-45713 _ 99 
; 5.3943 





4, = 0433022 = 05110 
4.5713 


O,, = tan" ae 0.0° 
; 0.2041 


8, = om (Se 





= 68.2250 
0.2041 


_ 0.2041 
cos0O 


= 0.2041 





3t 


= OAT 25502 


3 60868.2250 


2041 
ac ocieee 


1.1718 


0.5502 
Yys, === = 0.5083 


1.1718 





Rotor Blade Geometry 


H, =5.4055 — 3.5063 = 1.8992 in 





Ge Oe 24 3e07 ih 
Z, = PTCA IID 59:6 BSBA blades 
1.5827 


Choose Z» to match transonic compressor rotor. Therefore, Z. = 22. 
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1.8992 


Rev = Dg 4559-152 02818 


AR 








Stator Blade Geometry 


H, =5.3943 —3.7483 =1.6460 in 


C= 1.6460 


3 


= 1.3717 in 





._ 2-m-4,5713 


, -1.25 = 26.1740 blades 
, 1.3717 


Choose Z;3 to match transonic compressor rotor. Therefore, Z3 = 27. 





1.6460 
ARnos = 9g 4,5713- 105 7 12379 
27 
C= bOH00 5907 in 
“1.2379 
= 1.3297 = 1.0638 in 





3 
The design calculations have to iterated until the losses converge to a chosen criteria. A criterion 


of 0.0001 was chosen for both the hand and computer calculations. 
5s Iteration #2 


Only the results will be shown for all subsequent iterations. 


ROTOR STATOR 
Boo | 56.560435987322300/Bo 27.347157421794000 


B 
C 0.58601093209171 Ca 1.470018487969480 


L 
Cp; 0.01 4452454402386|Co, | 0.071059015271782 


0.06123341323541 
1.97841 707024513 


-0.02634785052288 1 


0.025691775310535 0.02546133547691 
0.05421 1496886609|or, | : 
0.957618541339988 : 


1.17184348697251 
1.13018217397965 





—| © — 
| Oo (o) 
©} 0 Oo 
| D oO 
oO] ® ~N 
N] © © 
NY] A © 
ol N oO 
o; = ~N 
Ny] © ie) 
| N N 
ao; — (o) 
Co} MO © 
Ol] & ice} 
OO} & £& 
oO] OO] CO] B® O}O| © 
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1.66081268021 1420/P/Pu | 1.468607319851730 
1.1781430790241 50am | 0.17223792203128 
1.774978623835910 0.696539265527826 
1.731530543151 750/As (in) | 48.558787505362500 
1.6993415116051 50|tam (in) | 4.57128697268207 
1.091667880043820|h; 


(in) 
(in) 
[saeraserzzareseny 
5597668878861. 
[eseaesmasorasay 
/aseaasesroosreaty 
/sara0eTa56900880 

ae ea 


3.493738020586840 
Table C1. Stage Performance Results (Iteration #2) 





0.3257971 82731520) Xo 


0.209326572089755] Xwor 





Table C2. Rotor Results (Iteration #2) 


TIP 


0.000000000000000)xs | 0.000000000000000 
0.352931481111 029/Xusr 0.513071730606907 
0.000000000000000\as = 0.000000000000000 
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0.429316175796641 0.429316175796641 


0.9091 00497709792|/Mwst | 1.326087214551410 
0.429316175796641 0.429316175796641 


Table C3. Stator Results (Iteration #2) 





ROTOR STATOR 


| CCSCROTORTs—“‘“OU*sLSOCOCtC(NNSTATOR 
1.924283336313540 1.690633536728720 
1 -603569446927950/Cs, | 1.408861280607270 
26.538058032755800|Z,' | 25.483562094627500 


22.000000000000000/Z, | 27.000000000000000 
0.994797734160298/AR' 1.271407814954980 


1.934346320096740 1.329733478780440 
1.272596263221 540/S3 | 1.063786783024350 


Table C4. Blade Geometry Results (Iteration #2) 





Compare , (iteration #2) and @, (iteration #1). 





O.,  — Or | = 0.05421 —0.02569 =0.02852 





Oy On. | = 0.08669 —0.02546 = 0.06123 
6. Iteration #3 
ROTOR STATOR 


56.560435987322300|Bo | 27.34715742179400 
Cc 


0.586010932091710}C. 1.47001848796948 


0.014264158667470 ; 


0.027683925940807|asetcs 0.059618090727575 
0.90877067413360 1.97841707024513 
0.99974659687064 1.00578807270894 


0.00047035065516 : 


0.025691775310535 0.02546133547691 
0.05384605190650 0.085079426204491 


0.95761854133998 i 


0.859357579736388|Ti/Ti | 1.17184348697251 


1.155975840871 94/Ts/T : 


OSFTC2 


ie) 


Pre/Pii 


| 


oa 


Co] N 


1.660812680211420|P3/P: 
1.1781430790241 50am 0.17223792203128 


1.774978623835910 : 
1.731823430894990 : 
) 


in 
1.699628954585760|ram (in) | 4.57128697268207 
1.091667880043820 : 

in) 
in 


PRAT TP ETT 


G 


0.219098107941022/ra ( 5.41615515589086 
0.772656848887 1 18/rs), (in) 3.72641 878947328 
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1.3262559555691 00/rnis 0.688019209608545 





[~snns86r080727675 
[o.nes0rer3e204a7 
[—v.seeo 9208800545 


53.865275938523100 
4.455879688743610 
1.923957899554990 


0.64488591 7495964 
5.417858638521 100 


3.4939007389661 10 
Table C5. Stage Performance Results (Iteration #3) 





[wage 22EEES 
Yon [| __O.300955580009557|¥a 
0.367525441 306733 


Table C6. Rotor Results (Iteration #3) 





TIP 


|< Se | | 
Ka 


Xweh 0.4077400371 67485|/Xwa | 0.552141570422119 


0.188552229153336 0.188552229153336 
0.3766591 76775861 0.510053392901174 
0.429316175796641 0.429316175796641 
0.909196076420228 1.325959056752780 
0.429316175796641 0.429316175796641 


Table C7. Stator Results (Iteration #3) 





ROTOR STATOR 
1 .923957899554990/H3 1.689736366417590 





1.6032982496291 60|C3) 1.408113638681320 
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26.542546935338600|Z;' | 25.497092664122100 


22.000000000000000}Z, 27.000000000000000 


2 
0.994629492953864/AR' 1.2707331 15606990 


1 -934346320096740/C3 | 1.329733478780440 
1.272596263221 540/S3 1.063786783024350 


Table C8. Blade Geometry Results (Iteration #3) 


Compare @,, (iteration #3) and @, (iteration #2). 











Op — Dro, | = 0.05385 —0.05421 = 0.0003654 
Oy = O..| = 0.08508 — 0.08669 = 0.001615 
7. Iteration #4 


ROTOR STATOR 
56.560435987322300[Bo 27.347157421794000 
0.58601093209171 Ca 1.470018487969480 
0.014266571444569 0.06922 1229482487 


0.027688608666518 0.059649745121757 
0.908770674133602|y* 1.978417070245130 


so 0.000470350655165 
or 
O.059850734692278[o%. | _—0.08511 1080598674 
PeilPu 
etn) | S3.86sse2668066700[ | 
fan Or) | 4a6587O6BETASBTOL | 
he (in) : 0 

(in) 

in) 


re) | 
a 
ra [__Sai7e6nr2sc07e00 
fn fay_[_Sa8BBESATEREO| 


Table C9. Stage Performance Results (Iteration #4) 


eo 
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TIP 
Noch 0.259622158556950 0.167426509603247 


[9.5296 7002TEEBOO|G, 
[_ORS7a7ETOTAO6AO)xe [0 259800967637099 
Por 


Ba 
: Xot : 
[820555227 E9EAS Ki 


Table C10. —_ Rotor Results (Iteration #4) 





TIP 


PP 


Table C11. Stator Results (Iteration #4) 





22.000000000000000/Z, 27.000000000000000 
0.994631648397049 1.270744874208740 


1 -934346320096740/C3 1.329733478780440 
1.272596263221 540/S3 1.063786783024350 


Table C12. | Blade Geometry Results (Iteration #4) 





Compare @,, (iteration #4) and @,, (iteration #3). 





ee O..| = 0.05385 — 0.05385 = 0.000004683 
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@,,_ —@,___ |= 0.08511—0.08508 = 0.00003 165 


Both the rotor and stator losses satisfy the criterion after the fourth iteration. 


P 


TI, =—* =1.6678 
Py 
T, 
t.=—2=1.1718 
Ty 
i4l 
_ 1.6678 — 


; | 0.9158 or 91.58% 
1.17181 


C2. CODE RESULTS 


1h 


oi 


h 


Eid 
| 


a 


th 


| 


poi / = Ose seen aeeAoae 


= 
(er) 


= 
o 


> 


i ed 
fa s097e6Ba08z6A00[ | 
jo | PaBeasaserarsse0009 [| 
indat [7-o8000000000000q) [| 
<=. 
a! 
| —== 
<< 


>) 


Pi 14.690000000000000 
Te 517.800000000000000 
Ro | 53.350000000000000 


1.400000000000000 
Table C13. Inlet Results 





[| CMEANS TTP 

19.153113641 505600|Bom | 47.8081003986021 OO/Ba | 60.32272154334440 
52.844869054566300|o2m | 45.871861 712689900 on | 40.226189216790300I 
0.326016514230627 0.282758297008589| X> 0.257833958315095 


0.328294695162428| Xijom 0.420281736792325] Xuyot 0.512268778422222 
3.455094238614300 4.4231997311 66900|ra | 5.391305223719500 
0.259909266378080| Xoo 0.203022938914057| Xoo 0.166566531027495 
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0.539065878334816|Dom | 0.455000000000000|D>; 0.390436898079840 


0.604152503515955|tsiam | 0.758468045193244)raa | 0.837422718264700 
[fe rrr0essec67s00601 CS” 


Fis isto ar298r90q 
[Foon 520000000000000) 
[[Fres ar8s300000000000) 
[Rs 6870000 0000000) | 
[ [fs ss 000000000000) 


Table C14. Rotor Results 





TIP 


Po CCKUB | CMMEANT TTP 
z 
Pf «| 0000000000000 | 
[ne asaaraeareaaeri 
an T250000000000009) | 
fax | e-800000000000009 J 
Fee [ 1025800000000009 J 
Re | T0e7vo0000000009 | 


Table C15. Stator Results 
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ROTOR STATOR 


FTC2 0.01509250281617 
0.00041528649613 
0.02562191983027 
0.041 12970914258 


0.95761854133998 1.67470350264381 
0.859357579736388] Tr3/Ti1 1.170653666750960 


xe} 
i 


© 


are 


— 
0 


4 


| 


nani 1.15465694167410IT3/Ti, 1.12899235375811 
1.47522045859872 


=] 


ai/Pr 1.65419002608182 


; 1.17663673828117 0.172309288336835 
EV/Pu | 1.767048277556160|A3/A, 0.692686802248097 


E2/Pu 1.734357050220880)A3 (in) | 48.440887367635200 
1.70368617510983 

TH 1.090701 41222377 
1.33002761786920 
0.218975111752094 
1 0.769783662636916 
> (in) | 53.832415426741300 
Table C16. Stage Performance Results 


Oo 


=] 


a 
ca 
a 
i 


al 
oO 


J 
mR 
2 
| 


ma 
is) 


Oo 


NO 
~ 


>] Oo] v 
ran 


> 





ROTOR 
1.93621098510521 OH; == | =: 1.698305724049770} 
1.613509154254340 

26.1916791 06060200/2;' |= 25.192519448425300 

Zs 


22.000000000000000/Z, | 27.000000000000000 
1.007953705186150 1.27717753305525 


i 
ye) 


1.920932454678200|C, 1.32051234513877 
1.263771351761 970/S3 | 1.056409876111010 


3.02727361 2963370 -14.428782843997100 


6.757152202149710|i*am | = -1.724296488710650} 
10.291 562745234000|i*a == | =: 8.241975399750170} 
-42.100896728952500 
-24.52312306331370 
-19.98318800786580 


3 


ili 
N No 


* 


ie) 
= 


ne} 
3 


S 


-5.757137628140870! 5*3, 21.988386237067600 
-2.8331 875653466 10! 5*3, 8.68868582584208 


-1.152333894067850! 83 4.39171197908887 
Table C17. | Blade Geometry Results 
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Ratio/Efficiency 


APPENDIX D. COMPARISON OF RESULTS 
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45 
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Hand Calc. Results Code Results 9 Transonic Compressor 


Figure D1. =Annulus Comparison 
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Figure D2. Stage Performance Comparison 
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Inches 














2.50 
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1.50 
1.00 
0.50 
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s2h s2m s2t s3h s3m s3t 
Solidity 


Hand Calc. Results Code Results 5 Transonic Compressor 


Figure D3. Solidity Comparison 
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1.50 ' 
i - ' 
: i 7 ' 
0.00 
c2 C3 H2 
Chord/Height 
Hand Calc. Results Code Results O Transonic Compressor 


Figure D4. Chord and Blade Height Comparison 
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VZ2h VZ2m VZ2t V2h V2m Vat 
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Figure DS. Rotor Velocity Comparison 
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Figure D6. —- Rotor Flow Angle Comparison 
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rst2h rst2m rst2t D2h D2m Dat 
Deg. of Reaction/Diffusion 
STransonic Compressor 


Hand Calc. Results 5 Code Results 


Rotor Degree of Reaction and Diffusion Comparison 
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VZ3m VZ3t V3h Vv3m V3t 


Velocity 





Hand Calc. Results Code Results OTransonic Compressor 


Figure D8. Stator Velocity Comparison 
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